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1.0  INTRODUCTION 

This  investigation  of  multiple  ducted  streams  with  a  periodic  driver 
flow  was  conducted  to  determine  if  the  periodic  driver  flow  was  an 
effective  mechanism  to  improve  the  momentum  and  energy  exchange  and  to 
enhance  mixing  between  compressible  high  speed  ducted  streams.  The 
ducted  system  investigated  was  a  two-dimensional  planar,  constant  area, 
air-to-air,  supersonic-subsonic  ejector.  A  large  scale  fluidic  oscilla¬ 
tor  was  developed  to  generate  the  periodic  pulses  in  the  primary  driver 
flow.  A  theoretical  analysis  based  on  the  assumptions  of  a  one¬ 
dimensional,  quasi-steady,  inviscid  flow  was  developed  and  compared  with 
experimental  results.  An  extensive  experimental  program  examining  both 
steady  and  periodic  driver  flows  in  the  ejector  was  carried  out. 

1.1  STATEMENT  OF  THE  PROBLEM 

The  impetus  of  this  investigation  is  rooted  in  fluid-dynamic  prob¬ 
lems  encountered  in  the  development  of  high  energy  chemical  laser  systems. 
The  high  Mach  number,  low  pressure  streams  entering  the  laser  cavity  do 
not  readily  mix.  This  results  in  a  short  lasing  zone  and  restricted 
laser  power  output.  The  introduction  of  an  unsteady  flow  into  this  flow 
could  potentially  enhance  the  mixing  between  these  multiple  high  speed 
compressible  laser  cavity  streams  by  causing  large  disturbances  and 
vorticity  in  these  flows. 

The  high  temperature,  low  pressure,  corrosive  laser  cavity  flow  must 
be  pumped,  after  exiting  the  cavity,  to  the  prevailing  ambient  pressure 
condition.  As  discussed  later,  ejector-diffuser  systems  are  well  suited 


2 


for  this  task,  despite  the  fact  that  the  ejector-diffuser  is  a  heavy, 
bulky,  and  inefficient  addition  to  the  laser  system.  Walker  [l]1  dis¬ 
cusses  topics  involved  in  laser-ejector  system  integration  such  as  packag¬ 
ing  multiple  ejectors,  diffuser  requirements,  and  staging  considerations. 
The  opinion  was  expressed  by  Walker,  et  al . ,  that  advances  in  ejector 
packaging  and  investigations  of  the  effects  of  laser  system  parameters 
on  ejector  operation  could  lead  to  improved  ejector  system  performances; 
however,  any  significant  improvement  in  ejector  performance  is  most 
likely  to  be  the  result  of  advances  in  ejector  technology. 

The  overall  goals  of  this  research  program  were  to  investigate  how  a 
periodic  driver  flow  influenced  ejector  performance  parameters  and  whether 
a  periodic  driver  flow  in  an  ejector  significantly  enhanced  mixing  and 
momentum  exchange  between  the  ducted,  high  speed,  compressible  primary 
and  secondary  streams. 

The  investigation  had  the  following  specific  objectives: 

1.  To  develop  an  understanding  of  the  mechanisms,  processes, 
and  interactions  involved  in  ducted  multiple  stream  flows 
with  a  periodic  driver; 

2.  To  develop  and  evaluate  a  theoretical  model  for  such  a 
system; 

3.  To  determine  if  a  periodic  driver  flow  demonstrates  the 
potential  to  make  ejector  systems  more  compact  and/or 
efficient; 


f Numbers  in  brackets  refer  to  entries  in  REFERENCES. 
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4.  To  develop  a  large-scale  fluidic  oscillator  with  no  moving 
parts,  frequency  adjustability,  and  small  losses; 

5.  To  add  to  the  limited  technology  and  data  base  of  such 
flows. 

1.2  GENERAL  BACKGROUND  AND  LITERATURE  SURVEY 

Ejectors  are  non-mechanical  devices  used  to  pump,  compress,  mix,  and 
convey  fluids.  Ejector  applications  include  high  vacuum  systems,  ram¬ 
jet  combustors,  thrust  augmentation  nozzles,  nuclear  reactor  core  flow 
recirculation  pumps,  and  chemical  laser  pumping  systems  to  name  just  a 
few.  Ejector  system  configurations  range  from  one  axisymmetric  primary 
driver  flow  nozzle  concentrically  directing  its  flow  into  a  constant  area 
cylindrical  mixing  tube  to  a  bank  of  hypermixing  nozzles  directing  skewed 
driver  flows  into  a  variable  area  rectangular  duct  [1,2,3]. 

The  basic  mechanisms  of  ejector  operation  are  viscous  interaction 
and  mixing  between  the  primary  driver  and  secondary  induced  flows.  Thrust 
augmentation  and  the  pumping  and  compression  of  a  secondary  flow  are  the 
result  of  the  momentum  and  energy  transfer  from  the  primary  to  secondary 
flow  within  the  mixing  duct  of  the  ejector  [4].  During  ejector  start-up 
the  momentum  of  the  primary  driver  flow  is  imparted  to  the  fluid  within 
the  mixing  duct;  this  causes  a  flow  towards  the  open  end  of  the  duct. 

The  resulting  pressure  drop  at  the  inlet  end  of  the  mixing  duct  induces  a 
secondary  flow  to  enter  the  duct.  The  primary  and  secondary  flows  inter¬ 
act  and  mix,  ideally  reaching  a  fully  uniform  state  at  the  exit  plane  of 
the  ejector  mixing  duct. 
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A  number  of  analytical  methods,  differing  in  levels  of  sophistication, 
have  been  applied  successfully  to  ejector  flows.  Finite  difference  tech¬ 
niques  were  used  by  Gilbert  and  Hill  [5]  to  model  a  two-dimensional 
planar,  variable  area,  slot  nozzle  ejector.  Hasinger  [6]  has  applied  a 
one-dimensional  analysis  to  a  number  of  ejector  mixing  modes  including 
constant  area,  constant  pressure,  and  a  defined  pressure  profile  along  the 
mixing  duct.  Tai  [7]  has  developed  a  computer  code  that  for  given 
entrance  conditions  and  a  wall  static  pressure  distribution  determines  the 
ejector  wall  contour  by  finite  difference  techniques .  A  one-dimensional 
inviscid  analysis  was  applied  by  Addy  and  Mikkelsen  [8]  to  the  supersonic- 
subsonic,  constant  area  ejector.  This  analysis  is  based  on  the  method  of 
Fabri  and  Paulon  [9].  The  quasi-steady  flow  ejector  analysis  developed 
during  this  investigation  is  based  on  the  steady  flow  ejector  analysis  of 
Addy  and  Mikkelsen.  More  detailed  analyses  by  Chow  and  Addy  [10,11]  and 
Chow  and  Yeh  [12]  of  supersonic-subsonic  ejector  flow  use  the  method  of 
characteristics  to  determine  the  initial  supersonic  primary  flowfield  and 
a  one-dimensional  model  of  the  secondary  flow.  Viscous  effects,  variable 
area  inlet  sections,  and  special  consideration  of  low  secondary  mass  flow- 
rates  are  incorporated  into  these  analyses. 

Supersonic-supersonic,  constant  area  ejector-diffuser  flow  was 
modeled  with  a  one-dimensional  control  volume  approach  by  Mikkelsen, 
Sandberg,  and  Addy  [13],  The  one-dimensional  theory  predicted  maximum 
ejector  compression  ratios  15  to  20  percent  higher  than  the  corre¬ 
sponding  experimental  values.  The  supersonic-supersonic,  constant 
area  ejector  was  particularly  susceptible  to  secondary  flow  separation 
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which  requires  a  more  sophisticated  analysis  than  the  one-dimensional  theory 
used. 

One-dimensional  analyses  of  ejector  flow  are  limited,  in  a  strict 
sense,  to  constant  area  ducts.  Constant  pressure  and  specified  pressure 
distribution  mixing  modes  can  be  treated  by  one-dimensional  analyses  [6,8] 
but  rely  on  more  sophisticated  analyses  or  empiricism  to  obtain  the 
required  wall  contour.  Viscous  interaction  between  primary  and  secondary 
ejector  streams  can  noticeably  effect  the  secondary  mass  flowrate  [10,11,12] 
under  certain  flow  conditions  but  one-dimensional  analyses  do  not  consider 
or  predict  such  flow  field  details.  Despite  these  limitations,  theoreti¬ 
cal  predictions  based  on  the  simple  one-dimensional  approach  are  generally 
in  good  agreement  with  experimental  results. 

Extensile  bibliographies  of  ejector  literature  and  papers  have  been 
compiled  by  Addy  and  Mikkelsen  [8,13]. 

Ejector  pumping  systems  have  their  advantages  and  disadvantages  when 
compared  to  conventional  mechanical  systems.  A  detailed  comparison  of 
ejector  and  mechanical  pumping  systems  has  been  made  by  Baker  and  McDermit 
[14].  The  no-moving-parts  simplicity  and  the  simple  geometry  of  ejector 
systems  facilitates  protective  cooling,  lining,  or  coating  of  wetted 
ejector  surfaces.  Ejectors  can  be  fabricated  from  most  materials  without 
difficulty.  Thus  ejectors  often  see  service  where  high  temperature  corro¬ 
sive  fluids  are  involved,  such  as  in  chemical  lasers.  Ejectors  start 
quickly  and  respond  rapidly  to  adjust  to  changing  flow  conditions. 

Ejector  applications  which  require  large  secondary  flowrates  and  high 
pressure  recovery  are  accomplished  by  the  packaging  of  ejectors  in  parallel 
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and/or  series.  A  diffuser  is  generally  used  to  further  improve  the  pres¬ 
sure  recovery  of  an  ejector  system.  Such  ejector-diffuser  systems  tend 
to  be  heavy  and  bulky  and  to  require  considerable  quantities  of  high 
pressure  primary  flow  to  pump  the  secondary  fluid.  Optimum  performance 
of  ejector  systems  is  usually  limited  to  a  small  range  of  operating  con¬ 
ditions  and  requires  total  system  integration  to  achieve. 

This  investigation  is  based  on  the  premise  that  a  periodic  driver 
flow  may  enhance  the  modes  of  flow  interaction  beyond  those  which  exist 
in  a  steady  flow  ejector.  Turbulent  mixing  and  entrainment,  the  basic 
mechanisms  involved  in  steady  ejector  flows,  have  been  the  subject  of 
numerous  investigations  [15-22].  Experimental  evidence  indicates  that 
large  scale  eddies,  which  evolve  from  surface  waves  generated  by  turbu¬ 
lence  induced  instabilities  at  a  fluid  interface,  act  as  the  primary 
mechanism  of  entrainment  in  turbulent  shear  flows.  The  rate  of  entrain¬ 
ment  and  mixing  layer  growth  depend  directly  on  the  intensity  of  the  large 
scale  eddies  and  indirectly  on  factors  which  determine  the  frequency  of 
eddy  generation  and  the  speed  of  their  decay.  These  large  scale  struc¬ 
tures  in  the  turbulent  mixing  layers  have  often  been  found  to  be  coherent 
and  periodic  in  nature. 

Numerous  approaches  have  been  taken  to  enhance  mixing  by  exciting 
instabilities  in  the  shear  layer  of  air  jets.  Crow  and  Champaign  [23] 
used  loudspeakers  in  a  plenum  chamber  upstream  of  the  exit  of  a  round  jet 
discharging  into  the  atmosphere  in  an  attempt  to  force  latent  periodic 
structures  that  may  exist  in  the  turbulence  of  the  jet  mixing  layer. 
Increased  entrainment  rates  and  turbulence  levels,  particularly  at  the 
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periphery  of  the  jet,  were  observed.  The  virtual  origin  of  the  mixing 
region  was  shifted  approximately  two  jet  diameters  upstream  with  a  2% 
fluctuation  of  the  exit  velocity.  Becker  and  Massoro  [24]  investigated  the 
instability  of  an  axisymmetric  jet  with  external  acoustic  stimulation. 

Pure  tones  and  noise  produced  both  large  increases  in  the  mixing  layer 
spread  rate  and  large  ring  vortices  that  disintegrated  into  small  scale 
turbulent  eddies.  Glass  [25]  located  a  reflecting  surface  near  the  exit 
of  an  overexpanded  supersonic  jet.  Disturbances  which  were  believed  to  be 
generated  in  the  shock  cells  of  the  supersonic  flow,  were  directed  towards 
the  jet  exit.  Periodic  vortex  shedding  was  observed  at  frequencies  as 
high  as  10  kHz.  The  spread  rate  and  velocity  decay  of  the  jet  were  notice¬ 
ably  increased.  McCormack,  Cochran,  and  Crane  [26]  also  obtained  increased 
mixing  rates  and  a  shortened  potential  core  by  mechanically  vibrating  a 
two-dimensional  slot  nozzle  that  was  discharging  air  into  the  atmosphere. 

Fully  and  partially  pulsed  jets  were  used  in  a  number  of  investiga¬ 
tions  to  study  the  effects  unsteady  flow  had  on  jet  mixing  and  turbulence. 
Binder  and  Favre-Marient  [27]  used  a  rotating  butterfly  valve  to  pulse  an 
axisyronetric  air  jet;  RMS  velocity  deviations  from  the  mean  of  up  to  30% 
were  produced  at  frequencies  up  to  200  Hz.  Centerline  velocity  decay  and 
spread  rates  reached  their  asymptotic  values  much  faster  than  for  the 
steady  jet  and  the  intensity  of  the  turbulence  increased  more  rapidly  and 
remained  at  higher  levels  well  after  pulsations  had  died  away  downstream. 
Bremhorst  and  Harch  [28]  fully  pulsed  an  axisymmetric  subsonic  air  jet 
using  mechanically  driven  rollers  that  only  allowed  air  flow  during 
approximately  1/3  of  each  revolution.  The  entrainment  capabilities  of  the 


jet  were  significantly  increased.  Olivari  [29]  examined  an  axisymmetric 
jet  that  was  almost  fully  pulsed.  A  fluidic  back  pressure  amplifier  that 
used  a  rotating  cam  to  modulate  the  control  port  was  used  to  pulse  the  flow. 
Turbulence  intensity  increased  significantly  but  mixing  between  the  jet  and 
the  surrounding  fluid  was  not  significantly  increased  over  that  of  the 
steady  jet.  It  appears  that  the  Strouhal  number,  St  =  uid/V ,  is  a  critical 
parameter  in  the  above  studies.  Olivari,  et  al . ,  had  St  values  of 
8.25  x  10  3  and  22  x  10  3  but  the  other  studies,  cited  above,  had  St  values 
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from  69  x  10  to  800  x  10  .  In  general  an  increase  in  the  Strouhal  num¬ 

ber  causes  an  increase  in  the  entrainment  of  a  jet  and  a  more  rapid  decay 
of  the  jet  core  velocity. 

Viets  [30]  developed  an  alternative  approach  to  stimulate  jet 
entrainment.  A  fluidic  nozzle  was  used  to  induce  a  jet  to  oscillate  trans¬ 
verse  to  the  flow  direction.  The  jet  spread  rapidly  and  takes  on  a 
serpentine  appearance  prior  to  jet  breakup.  Platzer,  Simmons,  and  Bremhorst 
[31]  compared  the  effects  on  the  entrainment  of  the  surrounding  fluid  of  a 
transverse  oscillating  jet,  fully  pulsed  jets,  and  co-axial  jets  with  a 
partially  pulsed  core.  The  fully  pulsed  jets  entrained  the  most  ambient 
fluid  and  the  partially  pulsed  and  transverse  oscillating  jet  entrained 
approximately  the  same  amount.  All  three  unsteady  jets  entrained  substan¬ 
tially  more  flow  than  a  corresponding  steady  jet. 

Ejectors  have  been  the  subject  of  a  number  of  pulsed  flow  investiga¬ 
tions  [32-35].  The  primary  driver  flow  was  fully  pulsed  by  a  rotating 
valve  mechanism  in  each  of  these  studies.  "Slugs"  of  primary  flow  are  dis¬ 
charged  into  a  duct  producing  compression  waves  in  front  of  the  "slug"  and 
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rarefaction  waves  behind  it.  The  piston-like  wave  action  of  the  primary 
flow  "slugs"  effectively  induced  a  secondary  flow  into  the  duct.  The  normal 
mode  of  energy  and  momentum  exchange  in  these  "wave  energy  exchangers"  is 
ideally  isentropic  and  more  efficient  than  the  irreversible  mixing 
mechanism  operating  in  steady  flow  ejectors.  The  secondary-to-primary 
mass  flowrate  ratio  of  the  pulsed  flow  ejector  was  much  larger  than  the 
mass  flowrate  ratio  of  the  corresponding  steady  flow  ejector  in  certain 
cases.  Since  the  secondary  or  induced  flow  stagnation  pressure  and  the 
exit  plane  static  pressure  were  atmospheric,  nothing  is  known  of  the  pres¬ 
sure  recovery  capabilities  or  the  effects  of  the  back  pressure  on  the 
performance  of  this  pulsed  ejector. 

In  some  cases  steady  flow  ejectors  have  been  found  to  resonate 
under  certain  pressure  conditions.  Quinn  [36]  has  examined  the  effects  of 
these  screech  tones  on  ejector  performance.  It  appears  that  resonance  may 
have  fixed  the  frequency  of  shear  layer  vortices,  thus  strengthening  them. 
Increased  entrainment  and  mixing  between  the  secondary  and  primary  flows 
were  observed. 

Experimental  evidence  strongly  suggests  mixing  between  multiple 
ducted  streams  and  ejector  performance  can  be  improved  by  an  unsteady 
periodic  driver  flow.  All  but  a  few  of  the  investigations  described  above 
were  concerned  with  a  significantly  different  flow  situation  than  is  con¬ 
sidered  in  this  investigation;  however,  turbulent  mixing  and/or  entrainment 
are  essential  to  each  one. 
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2.0  THEORETICAL  ANALYSIS  OF  THE  CONSTANT  AREA  SUPERSONIC- 
SUBSONIC  EJECTOR  WITH  A  PERIODIC  DRIVER  FLOW 


A  one-dimensional,  quasi-steady  analysis  was  used  to  model  the 
unsteady  ejector  flow.  A  computer  code,  UCAE,  was  written  to  perform 
this  analysis  on  either  experimentally  determined  input  conditions  or  on 
a  default  set  of  conditions  including  a  periodic  sinusoidal  driver  flow. 

Figure  2.0-1  indicates  the  location  of  the  stations  and  control  volume 
used  in  the  overall  analysis.  The  analysis  assumes  that  uniform  super¬ 
sonic  primary  and  subsonic  or  sonic  secondary  streams  enter  the  mixing 
duct  at  Station  1,  interact  and  mix,  and  exit  the  duct  as  a  uniformly 
mixed  subsonic  flow  at  Station  3.  Part  of  the  quasi-steady  analysis  is 
based  on  the  method  of  Fabri  and  Paulon  [9]  for  the  steady  flow  analysis 
of  an  ejector. 

2.1  EJECTOR  FLOW  REGIMES 

Four  ejector  flow  regimes  are  considered  in  this  analysis.  These 
flow  regimes  are  characterized  by  the  secondary-to-primary  static  pres¬ 
sure  ratio  at  Station  1,  P  / P  , ,  and  whether  or  not  the  flow  conditions 

si  pi 

at  Station  1  are  independent  of  the  back  pressure.  The  "saturated  supersonic" 
regime,  SSR,  can  exist  when  Psl/Ppl  1  1-0  and  Ms  l  =  1.0.  If  Ptl/P  >  1.0, 
the  secondary  flow  expands  against  the  primary  as  these  streams  enter  the 
mixing  duct  and  the  minimum  secondary  flow  area  is  at  Station  1.  The 
"supersonic"  regime,  SR,  can  exist  when  P^/P  <  1.0  and,  as  discussed 
later,  when  the  secondary  flow  is  choked  at  an  aerodynamic  throat  formed 
between  the  primary  flow  and  the  duct  wall.  An  "approximate"  regime, 

APPR,  is  used  to  model  flows  which  are  ideally  in  the  supersonic  regime 
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but  due  to  computational  difficulties  are  most  practically  modeled  by 
assuming  P  / P  ~  1.0. 

The  three  regimes  mentioned  above  can  occur  only  if  the  back  pres¬ 
sure,  P  ,  is  sufficiently  low.  If  the  back  pressure  is  too  high,  the 
secondary  flow  cannot  be  choked  in  the  supersonic  and  saturated  super¬ 
sonic  regimes  and  the  Psl/P  =  1.0  condition  used  in  the  approximate 
regime  is  no  longer  a  good  assumption.  The  "mixed"  regime,  MR,  exists 
when  the  back  pressure  is  too  high  to  allow  the  other  regimes  to  operate. 
In  the  mixed  regime  the  subsonic  secondary  flow  conditions  at  Station  1 
are  dependent  on  the  back  pressure  imposed  at  the  exit  plane,  Station  3. 
SR,  SSR,  and  APPR  flow  conditions  at  Station  1  are  always  independent  of 
the  back  pressure. 

The  flow  characteristics  which  distinguish  one  flow  regime  from 
another  are  "break-off"  conditions.  By  changing  ejector  operating  condi¬ 
tions,  ejector  flow  can  "break-off"  from  one  flow  regime  to  another.  For 
example,  SSR  flow  breaks-off  with  MR  flow  at  M  =1.0  and  P  ,/P  >1.0, 

si  si  pi 

th e  "break-off"  conditions.  If  M  ,  is  somehow  decreased  from  its 
"break-off"  value  of  1.0,  the  flow  breaks-off  from  SSR  to  MR  flow.  The 
theoretical  analysis  determines  ejector  flow  parameters  at  break-off 
points  between  the  mixed  regime  and  back  pressure  independent  regimes. 

For  the  given  input  data,  the  solution  procedure  used  to  determine 
ejector  flow  conditions  at  Stations  1  and  3  depends  on  the  basic  charac¬ 
ter  of  the  flow  that  is  possible  at  Station  1.  Each  of  the  four  flow 
regimes  considered  in  this  analysis  requires  a  different  computational 
approach.  The  determination  of  the  flow  regimes  which  can  exist  within 
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the  ejector  is  the  initial  step  taken  in  the  analysis.  This  initial 
determination  ic^ires  consideration  of  the  possible  secondary-to-primary 
static  pressure  ratio  at  Station  1,  Psl/Ppl-  The  back  pressure  imposed 
at  the  exit  plane  of  the  ejector  is  also  used  to  determine  which  one  of 
two  possible  regimes  exist,  MR  or  a  back  pressure  independent  regime. 

If  Psl/P  >_  1.0,  the  regime  of  operation  is  either  the  saturated 
supersonic  or  mixed  regime.  In  either  case,  the  secondary  flow  expands 
against  the  primary  flow  as  both  streams  enter  the  mixing  duct.  The 
minimum  secondary  flow  area  is  at  Station  1  so  0  <  M  <1.  When  the 
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secondary  flow  is  choked,  M  =  1,  the  regime  is  saturated  supersonic; 
otherwise,  the  mixed  regime  must  exist,  M  <  1.  The  mixed  regime  occurs 
when  the  back  pressure  is  too  high  to  allow  secondary  choking. 

Both  the  supersonic  and  mixed  regimes  are  possible  when  Psl/Ppl  <  1-0. 
The  primary  flow  expands  against  the  secondary  as  each  stream  enters  the 
duct.  It  is  assumed,  after  Fabri  and  Paulon  [9],  that  the  primary  and 
secondary  streams  remain  distinct  during  their  initial  interactions  and 
that  if  the  back  pressure  is  sufficiently  low,  the  secondary  flow  is  choked 
at  an  aerodynamic  throat  formed  between  the  primary  flow  and  the  duct  wall. 
If  secondary  choking  does  occur  at  some  location.  Station  2,  within  the 
duct,  see  Fig.  2.1-1,  the  ejector  operates  in  the  supersonic  regime.  An 
ejector  with  an  unchoked  secondary  flow  is  operating  in  the  mixed  regime. 

SR,  SSR,  and  MR  flows  do  not  include  all  the  possible  input  conditions 
which  require  a  separate  computational  solution  procedure.  The  supersonic 
and  saturated  supersonic  regimes  "break-off"  into  the  mixed  regime  if  the 
secondary  flow  becomes  unchoked,  however,  the  saturated  supersonic  regime 


computationally  does  not  necessarily  “break-off"  into  the  supersonic  regime 
at  the  conditions:  Ms x  =  1  and  Ptl/P  =1.  If  a  supersonic  regime  flow 

is  considered,  one  would  expect  that  as  P  ,/P  ,  is  increased  the  location 
of  the  choked  secondary  flow,  Station  2,  would  approach  Station  1  as 
Psl/P  approached  a  value  near  1.  Thus,  Ms 1  should  approach  Ms 2  or  a 
value  of  1.  The  computational  procedure  for  the  supersonic  regime  solution 
involves  a  search  algorithm  in  which  M  is  successively  incremented  and 
used  as  a  calculation  parameter.  Typically  for  a  certain  set  of  input 
conditions,  a  value  of  P  ,/P  ,  =  1  is  encountered  with  M  <  1.  A 
supersonic  regime  solution  is  no  longer  obtained  and  the  flow  conditions 
are  not  at  supersonic-saturated  supersonic  break-off  values.  The  lowest 
value  of  M  at  which  this  occurs  is  defined  as  the  critical  secondary 
Mach  number,  M  t  ,  and  marks  the  break-off  of  another  computationally 
distinct  regime  within  the  supersonic  regime.  It  has  been  found  that 
supersonic  regime  solutions  may  exist  for  M  <  Ms i  <  1*  These  solutions 
all  have  values  of  P  ,/P  ,  very  near  to  1.0  which  is  also  the  break-off 

sl  pi 

condition  between  the  regimes  of  supersonic  and  saturated  supersonic  flow. 
It  appears  tnat  solutions  in  this  regime  can  be  well  approximated  by  assum¬ 
ing  P^/P  =  1.  This  "approximate"  regime  appears  to  be  the  result  of  the 
sensitivity  of  the  flow  at  these  conditions  and  the  need  for  greater 
numerical  accuracy.  The  computational  procedure  in  the  supersonic  regime 
involves  a  search  technique  and  requires  iterations  to  determine  a  solution 
by  the  Fabri  criterion  at  each  step  of  the  overall  search.  It  appears  that 
a  substantial  increase  in  the  overall  number  of  iterations  would  be 
required  to  obtain  an  accurate  SR  solution  that,  at  these  flow  conditions, 
probably  does  not  significantly  differ  from  the  APPR  solution. 
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I 

^  Theoretically  the  supersonic,  saturated  supersonic,  and  approximate 

j  regimes  all  break-off  into  the  mixed  regime  flow  when  the  imposed  back 

pressure  is  greater  than  the  theoretical  value  at  the  break-off  condi¬ 
tions.  However,  as  previously  mentioned,  when  the  back  pressure  is  lower 
than  the  break-off  value,  the  flow  conditions  at  Station  1  are  unaltered 
from  their  break-off  values  and  are  independent  of  the  back  pressure. 

Figure  2.1-2  is  a  sketch  to  scale  of  the  solution  surfaces  determined  by 
the  analysis  for  the  ejector  configuration  which  was  examined  experi¬ 
mentally.  Break-off  values  establish  the  loci  of  break-off  points  along 
which  the  SR,  APPR,  and  SSR  surfaces  intersect  the  MR  surface.  In  the 
analysis  and  in  Fig.  2.1-2,  Ps0/Pp0  is  an  input  flow  condition,  Pn)3/Pf0 
is  the  non-dimensional  back  pressure  boundary  condition,  and  M  is  the 
dependent  variable. 

2.2  OVERALL  CONTROL  VOLUME  ANALYSIS 

The  control  volume  used  in  the  overall  analysis  is  indicated  in 
Fig.  2.0-1.  It  was  assumed  that  piecewise  uniform  flows  enter  the  duct 
at  Station  1  and  exits  subsonically  at  Station  3  as  a  uniformly  mixed 
flow.  A  unique  solution  is  obtained  by  determining  from  the  input  flow 
conditions  which  of  the  flow  regimes  is  possible  when  the  back  pressure 
is  sufficiently  low.  Flow  conditions  at  Station  1  can  then  be  determined. 
If  the  back  pressure  is  greater  than  the  break-off  value,  the  mixed  regime 
exists  and  a  unique  solution  is  found  by  satisfying  the  back  pressure 
boundary  condition.  The  conservation  principles  are  applied  to  the  over¬ 
all  control  volume,  Fig.  2.0-1,  to  determine  both  the  exit  plane  conditions 

i 

I 


6 


at  the  break-off  point  and  the  mixed  regime  solutions.  The  overall  control 
volume  analysis  has  a  unique  solution  for  the  back  pressure  independent 
regimes,  i.e.,  the  break-off  solution.  The  analysis  cannot,  for  instance, 
find  the  exit  plane  Mach  number,  M  ,  when  the  back  pressure  is  less  than 
the  break-off  value. 

The  theoretical  analysis  with  an  unsteady  driver  flow  was  carried-out 
under  the  assumption  of  quasi-steady  flow.  The  periodic  input  flow  condi¬ 
tions  were  subdivided  into  time  increments  and  the  theoretical  analysis 
was  performed  during  each  time  increment  as  if  the  flow  were  steady  and 
independent  of  other  time  increments.  Approximate  integral  averages  of 
the  performance  parameters  can  then  be  obtained. 

The  following  assumptions  are  made  in  the  analysis: 

3/  ) 

1.  Quasi-steady  flow  3^'  =  ®  ^or  eac^  time  s^eP» 

2.  Piecewise  uniform  flow  at  Station  1; 

3.  Fully  mixed  uniform  flow  at  Station  3; 

4.  Frictional  effects  within  the  duct  are  negligible; 

5.  No  shear  or  shaft  work; 

6.  Gravitational  potential  energy  changes  are  negligible; 

7.  The  flow  is  adiabatic  and  non-reacting  within  the  duct; 

8.  Primary  and  secondary  gases  are  calorically  perfect; 

9.  The  primary  nozzle  base  area  is  negligible;  and 

10.  Flows  are  brought  isentropical ly  from  their  stagnation 

states  to  Station  1. 

Continuity  may  then  be  expressed  as: 

0  p  V  •  dA  =  0  ( 

‘'cs 


2.2- 
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p  V  A  +  p  V  A  =  p  _V  A  (2.2-2) 

pi  pi  pi  i  1  «1  >1  m3  it3  m3 

If  the  mass  flowrate,  w  =  p  AV,  is  employed,  the  continuity  equation 
may  be  written  as: 


w  +  w  -  w  . 

s  p  m3 


(2.2-3) 


w  +  1  =  w  ,/w 

s  p  m3  p 

where  w  h  w  /w  is  the  mass  flowrate  ratio. 

i  p  *  p 

The  mass  flowrate  function  is: 


w  R 


PA  Mw  o 


=  M  y  1  +  |- 


w 


5  fj (y»M) 


(2.2-4) 


(2.2-5) 


Expressing  the  mass  flowrate  ratio,  Wsp,  in  terms  of  the  mass  flowrate 
function  yields: 

P  ,  A  ,  [mw  T  11/2  f  fy  ,M  1 

_L  .  p°  1  1  »  » 1 )  /■>  ■> 


.p  '  p  .  A 

pi  pi 


Mw  T  „ 

p  >0 


fl 

[y  ,m  ,1 

«  J  1 

fl 

-j* 

(2.2-6) 


Conservation  of  momentum  in  the  X  or  flow  direction  may  be  expressed 


ZF*  "  0  V  (p  V  •  dA)  • 

;CJ  X 


(2.2-7) 


The  momentum  equation  for  the  overall  control  volume  is: 

P  ,A  ,  +  P  ,A  ,  -  P  ,A  =  p  TA  V2,  -  p  A  ,V2,  +  p  ,A  ,V2, 

plpl  s  1  J  1  m3  m3  m3  m3  m3  I  'p  1  p 1  p 1  ri  1  i  1  >  1 


(2.2-8) 


Using  the  relationship 
p  V2  =  PyM2 


(2.2-9) 
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Eq.  (2.2-8)  can  be  written  as: 

li!  .  Jli-  ( 1  +  if  I  *  fl  +  Y  «!  |  =  .  J!±  [l  +  y  H2  |  .  (2.2-10) 

pp,  l  '•  *>J  l  -)  P„  l  ■») 

Substituting  Eq.  (2.2-6)  for  P  /P  and  expressing  P  ,/P  ,  in  terms  of  the 

s  1  p  i  '  rr.3  pi 

mass  flowrate  function,  Eq.  (2.2-5),  produces: 


Mvmp. 

)  *  Mv 

Mw  T 

p  so 

Mw  *  T7 

s  pO 

1  /  2 

•  w 

s  p 

?|  % 
It- 

• 

3°4  2 

i  i 

1/2 

C +  ”.p) 

f,  (y  ’M  J 

2  m  m3 


where  the  function  f2(y,M)  is  defined  as: 


f,(y,M) 


1  +  yM 


M  y  1  +  -  M2 

Solving  Eq.  (2.2-12)  for  M  produces 


1/2  ‘ 


- 

f  j2 

K  - 2 

L 2 _ 

+ 

a 

f2  *  2 

2 

+  2 

[y-i] 

Y 

if2  2y 

r*  ‘  wrj. 

1/2 

<Y-l)p  -  j 

Jy 

Y-Dj 

Only  the  subsonic  solution  of  Eq.  (2.2-13)  is  used. 
The  energy  equation  can  be  written  in  the  form: 

nin  0W  r  [  .,2 

gZ  (pV  •  dfl)  . 

''CS  ' 


(2.2-11) 


(2.2-12) 


(2.2-13) 


(2.2-14) 


Under  the  foregoing  assumptions,  the  energy  equation  simplifies  to: 

ji  h0(p  V  •  dA)  =  0  (2.2-15) 

'Cl % 


where  h  =  h  +  V  / 2.  For  adiabatic  flow,  the  overall  control  volume 


analysis  energy  relation  is: 


(2.2-16) 


Dalton's  law  of  partial  pressures  is  assumed  to  be  applicable  to  the 
mixing  process  within  the  control  volume  as  a  means  of  determining  the 
mixed  flow  gas  properties  at  Station  3.  Mass  flowrate  ratios  are  used  to 
express  the  mass  fractions  in  the  relations  for  a  constant  volume  mixing 
process.  This  is  perhaps  the  most  unrealistic  part  of  the  quasi-steady 
approach  for  non-constant  property  flows.  The  properties  of  the  "slugs" 
of  flow  currently  entering  the  duct  are  used  to  determine  the  exit  plane 
properties  at  the  same  point  in  time.  These  conditions  at  Station  3  are 
then  assumed  to  be  the  back  pressure  conditions  affecting  the  flow  at 
Station  1.  If  amplitudes  and/or  phase  lags  are  small  this  will  not  result 
in  a  significant  error.  The  mixed  flow  properties  expressed  in  terms  of 


w  and  the  primary  and  secondary  gas  properties  are: 

S  p 


(2.2-18) 


(2.2-19) 
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Using  the  relationship 

r  -  yR 
p  '  Mw(y- 1 )  ■ 

the  ratio  of  specific  heats  can  be  expressed  as 
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(2.2-20) 


(2.2-21) 


Thus  the  mixed  flow  gas  properties  are  defined  in  terms  of  ys  ,  y  ,  Mws  , 
Mw  ,  and  w  . 

p  « p 

The  non-dimensional  form  of  the  back  pressure  boundary  condition, 
P„n3/Ps0  can  be  found  using  Eqs.  (2.2-4)  and  (2.2-5).  The  result  is 
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(2.2-22) 


The  preceding  relations  are  used  to  formulate  a  solution  procedure 
for  the  back  pressure  independent  and  mixed  flow  regimes.  In  this  pro¬ 
cedure,  it  is  assumed  that  fy  ,  y  ,  M  ,  ,  P  ,/P  .  ,  A  , /A  , ,  Mw  /Mw  , 

T  „/T  ace  known.  In  the  mixed  regime,  unknowns  include  (w  ,  P  ,  /P  ,  , 

sOpOj  3  lip’ll'  pi’ 

Mm3  ’  f  ii’  ^mo^po)  •  The  problem  for  the  back  pressure  independent  regimes 
has  some  of  these  variables  specified  from  the  outset. 
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In  the  saturated  supersonic  regime,  M  =  1.0;  in  the  approximate 
regime,  P  ,/P  ,  =1.0.  The  static  pressure  ratio,  P  /P  ,  can  be 

slpl  si  pi 

determined  for  a  value  of  M  by  the  Fabri  criterion  in  the  supersonic 
regime;  in  this  regime,  M  is  used  as  a  calculation  parameter.  As  pre¬ 
sented,  the  problem  is  not  closed.  There  are  five  mixed  regime  unknowns 
but  only  three  conservation  equations  and  the  back  pressure  boundary 
condition  are  available  to  determine  them.  The  back  pressure  independent 
regimes  have  four  unknowns  and  three  equations.  One  basic  unknown  must 
be  specified  as  a  calculation  parameter.  The  secondary-to-primary  stag¬ 
nation  pressure  ratio  has  been  cnosen  as  the  calculation  parameter.  The 
pressure  ratio,  Pj0/Pp0,  is  fundamental  to  both  flows  and  can  be  experi¬ 
mentally  measured  versus  time  without  difficulty.  Thus  P^/P  can  then 
be  determined  from  the  isentropic  relation: 


-i) . 


fy -11  ,  V“tp 

I1  *  hrk. 


(2.2-23) 


for  a  given  M  ;  the  problem  is  closed. 


2.3  SUPERSONIC  REGIME  ANALYSIS 

A  control  volume  analysis,  see  Fig.  2.1-1,  which  is  based  on  the 
method  of  Fabri  and  Paul  on  [9]  has  been  used  to  determine  supersonic 
regime  secondary  flow  conditions  at  Station  1.  Due  to  the  interaction 
between  the  primary  and  secondary  streams  within  the  duct,  the  assumption 
is  made  that  the  secondary  flow  chokes  within  the  duct  at  Station  2  and 
that  the  secondary  and  primary  flows  separately  co-exist  to  fill  the  duct. 
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The  pressure  ratio,  P^/P  is  found  for  an  assumed  value  of  Mf x  by  apply¬ 
ing  the  conservation  of  mass  and  momentum  between  Stations  1  and  2.  M  is 
determined  by  iteration  by  a  search  routine  until  the  value  of  P  /P 

1 1  p  i 

obtained  from  the  Fabri  criterion  agrees  with  the  value  resulting  from 
input  conditions:  M  , ,  P  „/P  and  the  assumed  M  , . 

pl  sO  pO  i  1 

The  following  assumptions  have  been  made: 

1.  The  primary  and  secondary  streams  remain  distinct  and 
uniform  between  Stations  1  and  2; 

2.  The  mean  pressures  of  the  primary  and  secondary  flow 
may  differ  at  any  cross  section; 

3.  P  ,/P  <  1.0; 

si  pi 

4.  M  =  1.0; 

s  2 

5.  The  flow  is  isentropic; 

6.  The  body  forces  are  negligible; 

7.  The  wall  friction  can  be  neglected; 

8.  No  shear  or  shaft  work; 

9.  The  primary  nozzle  base  area  is  negligible;  and 

10.  The  flows  are  brought  isentropically  from  their  stagna¬ 
tion  states  to  Station  1. 

The  area  ratio,  As1/Aj2,  can  be  expressed  with  the  area  ratio  function, 

A/A*, 

A/A*(y,M)  =  1 
as 

\  l  /A,  2  =  A/A*(ye  ,Ms,) 
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=  fj(y.M) 


(2.3-1) 
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since  Mj2  =  1.  The  flow  area  available  to  the  primary  at  Station  2  and 
the  determination  of  Mp2  follow  directly  from  the  conservation  of  mass; 
the  relationship  is 


or 
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(2.3-2) 


M  , ,  found  numerically,  must  be  greater  than  M  . 

P  2  pi 

Conservation  of  momentum  between  Stations  1  and  2  produces: 


P.,A.,(’  *  *."?,)  +  P,.Api('  ♦*,*,)  -  P2A2(1  +Y.) 

*  pPA,(' +  \K*} 


(2.3-3) 


Dividing  by  PplApl  and  rearranging  for  computational  simplicity,  the 
expression  for  P  , /P  ,  is: 
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(2.3-4) 


With  the  input  values  of  (y  ,  y  ,  M  ,  ,  A  ,/A  _)  and  an  assumed  trial 

l  *  p  pl  pi  m3  J 

value  of-M$1,  Eq.  (2.3-4)  can  be  solved  after  determining  M  from 

Eq.  (2.3-2)  and  the  isentropic  relations.  The  supersonic  regime  solution 
for  flow  conditions  at  Station  1  has  been  completed  when  a  value  of  M  i  is 
found  that  results  in  convergence  of  the  values  of  P  ,/P  ,  determined  from 

*  1  p  1 

Eqs .  (2.2-23)  and  (2.3-4). 


2.4  COMPUTATIONAL  PROCEDURE 


The  steps  in  the  computational  procedure  are: 

1.  To  determine  which  back  pressure  independent  regime  can 
exist, if  the  back  pressure  is  sufficiently  low,  at  the 
assumed  value  of  the  calculation  parameter,  Ps0/Pp0'» 

2.  To  use  the  appropriate  Station  1  conditions  and  the 
overall  control  volume  analysis  to  determine  the  result¬ 
ing  exit  plane  break-off  conditions;  and 

3.  To  perform  an  iteration  for  a  mixed  regime  solution  if 
the  back  pressure  condition  imposed  at  the  exit  plane 
is  greater  than  the  break-off  value  determined  in  item 
2  above. 


The  first  step  taken  to  determine  which  back  pressure  independent 
flow  regime  can  exist  is  to  determine  the  flow  conditions  at  which  the 
supersonic  and  approximate  flow  regimes  break-off  from  one  another.  Mcrjt 
the  break-off  value  between  these  regimes  of  the  secondary  flow  Mach 
number  at  Station  1,  is  found  by  carrying  out  the  supersonic  regime 


analysis  in  Section  2.3  with  assumed  values  of  M  .  The  minimum  value  of 

s  1 

M  ,  which  produces  P  ,  /P  ,  =1.0  from  Eq.  (2.3-4)  is  the  critical  Mach 

Si  r  Slpl  * 

number,  M  .  The  search  for  M  ,  is  not  difficult  because  P  ,/P  , 

crit  crit  slpl 

increases  monotonically  with  M  in  the  supersonic  regime.  The  correspond 


ing  break-off  value  of  Ps0/Pp0  is  found  using  the  relation 


P  ,/P  „ 

pi  P° 

P  ,/P  „ 

si  s  0 


(2.4-1) 


where  the  isentropic  relations  and  P  ,/P  ,  =1.0  are  used. 

slpl 
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The  possibility  of  a  SSR  flow  can  be  determined  by  assuming  M  x  =1.0 
and  solving  for  the  resulting  value  of  Ptl/P  .  The  relation  is 
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If  p5l/ppl  :L  1.0,  either  the  saturated  supersonic  or  mixed  flow  regime  must 
exist.  When  P  ,/P  ,  <1.0  and  p„/p„>fp„/p„)  ,  the  flow  regime 

must  be  either  approximate  or  mixed.  Either  the  supersonic  or  the  mixed 
regime  exist  for  values  of  P  „/P  „  lower  than  fP  „/P  1  .  Thus,  as 

s  o'  pO  I* O'  pO  J 

v  J  c  r  i  t 

previously  discussed,  back  pressure  independent  flow  conditions  at  Station 
1  are  determined. 

The  overall  control  volume  analysis  is  then  used  to  determine  the 
flow  conditions  at  Station  3.  This  is  accomplished  by  solving  in  order 
Eqs .  (2.2-23),  (2.2-6),  (2.2-20),  (2.2-21),  (2.2-18),  (2.2-19),  and  (2.2-17) 
for  P  ,/P  ,  ,  w  and  mixed  flow  gas  properties.  P  , /P  ,  =  1.0  is  assumed 
for  APPR  flow  so  Eq.  (2.2-23)  is  used  to  determine  M  .  M  .  is  found  by 
solving  Eqs.  (2.2-11)  and  (2.2-13).  Only  the  subsonic  solution  of 
Eq.  (2.2-13)  is  considered.  The  break-off  back  pressure  ratio,  P  ,/P  . 

is  determined  by  using  Eq.  (2.2-22).  If  the  imposed  back  pressure  ratio 
is  greater  than  this  break-off  value,  the  mixed  regime  exists.  In  the 
mixed  regime,  an  iterative  procedure  is  necessary  to  find  the  value  of  Mgl 
which  causes  P^j/Pl0  to  converge  to  the  given  value.  Each  iteration  for 
a  MR  solution  must  start  with  Eq.  (2.2-23)  and  then  proceed  with  the  over¬ 
all  control  volume  analysis  as  discussed  above. 

A  computer  code,  UCAE,  was  written  to  perform  this  analysis.  Section 
6.0,  the  APPENDIX,  contains  UCAE  and  a  sample  input  and  output. 


2.5  PARAMETRIC  STUDY 


A  parametric  study  of  theoretical  ejector  results  was  made  with  UCAE 
for  a  wide  range  of  P  /P  „  values.  The  input  conditions  used  for  this 

s  0  p  0  r 

parametric  study  correspond  to  the  experimental  setup.  Values  of  w  , 

S  p 

P  , /P  , ,  M  , ,  and  P  /P  „  at  the  break-off  conditions  of  the  back  pres- 
sure  independent  regimes  with  the  mixed  regime  are  plotted  versus  P.0/Pp0 
in  Figs.  2.5-1,  2.5-2,  2.5-3,  and  2.5-4.  The  lowest  value  of  Ps0/Pp0 
plotted  in  these  figures  corresponds  to  M  =  .001,  the  minimum  value 
considered  by  UCAE. 

The  mass  flowrate  ratio.  Fig.  2.5-1,  is  a  smooth  monotonic  function 
of  ps0/pp0>  increasing  with  Ps0/Pp0  in  all  three  back  pressure  independent 
flow  regimes.  The  slight  change  in  slope  between  the  SR  and  SSR  segments 
of  the  curve  occurs  gradually  over  the  APPR  segment  of  the  curve. 

The  static  pressure  ratio,  P  ,/P  , ,  is  a  linear  function  of  P  /P  „ 

r  si  pi  sOpO 

for  SSR  flow  but  is  nonlinear  for  SR  flow,  see  Fig.  2.5-2.  Discontinuities 
exist  in  the  P  ,/P  ,  solution  at  the  SR-APPR  break-off  and  the  APPR-SSR 
break-off.  The  SR-APPR  discontinuity  is  due  to  the  Ptl/P  =  1.0  approxi¬ 
mation  for  APPR  flow  and  insufficient  numerical  accuracy,  as  previously 
discussed.  The  APPR-SSR  discontinuity  is  due  to  the  constant  secondary 
choked  flow  area  in  SSR  flow  which  imposes  an  additional  constraint  on 
the  nature  of  the  solution.  This  is  more  apparent  in  the  plot  of  M 
versus  ps0/pp0>  Pig-  2.5-3. 

Figures  2.5-4  and  2.5-1  shows  the  relationship  existing  between 
ejector  pressure  recovery  and  mass  flowrate  performance.  Large  mass  flow- 
rates  are  obtained  in  the  saturated  supersonic  regime  at  the  sacrifice  of 


the  pressure  recovery.  The  converse  is  true  for  supersonic  regime  flows  at 


lower  values  of  P  „/P 

s  0  pO 

Figures  2.5-3  and  2.5-4  are  the  normal  projections  of  the  loci  of 
break-off  points  between  the  mixed  and  back  pressure  independent  flow 
regimes,  as  shown  in  Fig.  2.1-2,  onto  the  |Pj0/Pp0,  Mjl|  and 
(ps0/pPo’  P,n3/Pso)  p1anes’  respectively. 
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Figure  2.0-1  Control  volume  used  in  the  overall  analysis 


Theoretical  mass  flowrate  characteristics  of  the  back  pressure  independent  flow 
regimes  versus  P  „/P 


Figure  2.5-3  Theoretical  values  of  M  for  the  back  pressure  independent  flow  regimes 
versus  P  /P  „  ‘ 
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3.0  THE  EXPERIMENTAL  PROGRAM 

A  scries  of  small-scale  cold  air  flow  experiments  were  conducted  at 
the  University  of  Illinois  Mechanical  Engineering  Laboratory  to  determine 
what  effects  a  continuous  periodic  driver  flow  had  on  the  performance  and 
flow  characteristics  of  a  two-dimensional,  planar  supersonic-subsonic 
ejector.  A  wide  range  of  operating  conditions  were  examined  with  both 
steady  and  unsteady  primary  driver  flows. 

3.1  EXPERIMENTAL  APPARATUS 

Filtered  and  dried  compressed  air  at  a  maximum  pressure  level  of 

3 

850  kPa  was  supplied  from  a  compressor  connected  to  a  100  m  storage 
tank  farm.  Steady-state  operation  was  possible  for  experiments  requiring 
mass  flowrates  of  approximately  1  kg/s  or  less  of  compressed  air.  This 
air  supply  system  was  the  source  of  primary  flow  in  all  ejector 
experiments. 

The  continuous  flow  test  facility,  which  was  used  to  run  these 
experiments,  consists  of  two  test  stands  and  a  central  control  panel. 

Each  test  stand  has  a  two  and  a  three  inch  compressed  air  supply  line. 

The  stagnation  pressure  in  either  line  can  be  adjusted  independently  at 
the  control  panel  by  air  operated  control  valves.  Two  sets  of  U-tube 
manometers  on  the  control  panel  are  used  to  measure  pressure  differences 
across  VDI  standard  nozzles  to  determine  mass  flowrates.  Control  panel 
switching  allows  the  selection  of  U-tube  manometers  with  fluid  specific 
gravities  of  1.0,  2.95,  or  13.7.  A  Wallace  and  Tiernan  Series  1500 
0-150  psia  pressure  gage  with  a  .0667%  full  scale  accuracy  was  used  to 
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measure  pressures  at  various  line  and  test  locations  and  to  calibrate  some 
pressure  transducers.  A  Wallace  and  Tiernan  precision  manometer  and  a 
control  panel  regulated  vacuum  source  were  also  used  for  pressure  trans¬ 
ducer  calibration.  Each  test  stand  has  20  pressure  taps  connected  to  a 
40  tube  mercury  manometer  board.  The  manometer  board  was  used  to  measure 
the  secondary  stagnation  pressure  and  to  obtain  a  quick  visualization  of 
the  approximate  static  pressure  distribution  within  the  ejector  mixing 
duct. 

The  ejector  primary  flow  was  taken  from  a  2-inch  supply  line  on  one 
of  the  test  stands.  A  VDI  standard  nozzle,  located  upstream  of  the  test 
stand,  was  used  to  measure  primary  mass  flowrates.  The  supply  line  was 
reduced  to  1%  inch  pipe  at  the  test  stand  source.  A  hose  fitting  was 
machined  from  a  standard  1^  inch  pipe  nipple  so  that  nylon  reinforced 
flexible  tubing  could  be  used  to  connect  the  primary  flow  directly  to  the 
ejector  inlet  for  steady  experiments  or  to  the  inlet  of  a  large-scale 
fluidic  oscillator  for  the  unsteady  driver  flow  experiments.  Figure  3.1-1 
is  a  schematic  of  the  experimental  apparatus  which  details  primary  and 
secondary  flow  lines. 

The  secondary  flow  was  induced  through  a  smooth  elliptical  inlet 
nozzle  from  the  atmosphere.  A  VDI  standard  nozzle  block,  located  in  the 
middle  of  a  2.13  m  length  of  seamless  2-inch  tubing,  was  used  to  deter¬ 
mine  the  secondary  mass  flowrate.  The  flow  was  throttled  with  a  gate 
valve  to  the  desired  stagnation  pressure  after  it  had  passed  through  the 
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VDI  nozzle  and  a  .35  m  holding  tank.  A  four  port  manifold  was  used  to 
distribute  the  flow  to  the  secondary  inlets  on  the  ejector  unit.  The 


threaded  holes  seen  in  Figs.  3.1-2  and  3.1-3  are  the  secondary  flow  inlet 
ports. 

The  holding  tank,  mentioned  above,  was  used  to  avoid  mass  flowrate 
measurement  errors  typically  encountered  when  a  metering  orifice  or 
nozzle  is  used  in  the  measurement  of  an  unsteady  flow.  The  deviation  due 
to  unsteady  or  pulsatile  flow  of  the  metering  device  discharge  coeffi¬ 
cient  from  steady  flow  values  is  not  well  understood.  Investigations 
[37,38]  have  shown  that  frequency,  amplitude,  wave  form,  and  the  pressure 
difference  across  a  flow  metering  device  all  have  to  some  degree  an  effect 
on  the  discharge  coefficient  in  a  periodic  flow.  A  decoupling  of  the 
dynamic  pressure  difference  across  the  metering  device  with  the  velocity 
through  it  as  well  as  vortex  shedding  phenomena  due  to  the  unsteady  flow 
are  suspected  of  erroneously  altering  mass  flowrate  measurements. 

Pressure  drops  across  a  metering  nozzle  or  orifice  in  a  periodic  flow  are 
not  representatively  measured  by  U-tube  manometers. 

For  experiments  conducted  with  pulsating  flow,  the  holding  tank  adds 
sufficient  capacitance  to  the  secondary  flow  system  so  that  pulsations 
are  damped  in  the  flow  to  a  negligible  level  at  the  VDI  standard  nozzle. 
The  holding  tank  that  was  used  has  more  than  adequate  volume  to  damp  the 
flow  according  to  the  Hodgson  number  criterion  discussed  by  Ower  and 
Pankhurst  [39] . 

Pulsation  of  the  primary  flow  was  achieved  with  a  large-scale  fluidic 
oscillator.  No  moving  parts,  simplicity,  small  losses,  and  frequency  and 
amplitude  variability  were  reasons  for  using  a  fluidic  device  to  produce 
the  pulsations  in  the  driver  flow.  An  initial  effort  during  this  program 
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involved  the  selection  and  development  of  a  fluidic  oscillator  configura¬ 
tion  to  produce  the  pulsating  primary  ejector  flow.  Scale  models  of 
promising  fluidic  oscillator  geometries  were  fabricated  from  styrofoam 
for  performance  testing  using  the  hydraulic  analogy.  The  hydraulic 
analogy  is  a  convenient  way  to  study  and  visualize  unsteady  compressible 
flow;  thus  it  provided  a  qualitative  means  to  understand  the  factors 
influencing  the  oscillator  performance  [40].  The  particular  oscillator 
configuration  used  in  the  air  flow  experiments  was  chosen  because  of  its 
capability  to  self  excite,  its  pressure  level  insensitivity,  frequency 
stability,  its  large  amplitude  pulses,  and  the  ease  with  which  its  operat¬ 
ing  frequency  can  be  varied.  Essentially  the  same  oscillator  configuration 
was  used  by  Halbach,  Otsap,  and  Thomas  in  a  fluidic  temperature  sensor  [41]. 

This  oscillator  consists  of  a  fluid  jet  which  impinges  on  a  knife- 
edge,  a  feedback  loop,  and  a  resonance  chamber;  the  body  of  the  oscillator 
is  shown  in  Fig.  3.1-4.  A  jet  impinging  on  a  knife-edge  is  in  a  bi-stable 
condition;  the  jet  tends  to  flip  to  one  side  of  the  knife-edge  or  the 
other.  This  flipping  action  results  in  a  pressure  pulse  moving  through 
the  portion  of  the  oscillator  fed  by  the  jet.  Pulses  are  either  reflected 
back  to  the  jet  at  the  end  of  the  resonance  chamber  or  looped  back  to  the 
jet  through  the  feedback  loop;  in  either  case,  the  interaction  with  the  jet 
causes  the  jet  to  flip  to  the  other  side  of  the  knife-edge.  The  resonance 
chamber  has  an  opening  at  the  reflecting  end,  as  shown  in  Fig.  3.1-5, 
through  which  the  periodic  primary  driver  flow  enters  the  ejector.  Flex¬ 
ible  tubing  made  up  a  part  of  the  resonance  chamber  between  the  ejector 
inlet  and  the  oscillator  body.  Similarly,  the  feedback  loop  was  completed 
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outside  of  the  body  of  the  oscillator  with  flexible  tubing  and  a  machined 
180°  turn.  Oscillator  frequencies  were  adjusted  by  changing  the  lengths 
of  the  resonance  chamber  and  the  feedback  loop  tubing.  Stagnation  pres¬ 
sure  probes  were  located  in  the  hose  fitting  connecting  the  source  flow 
to  the  oscillator  inlet  and  in  the  hose  fitting  at  the  reflecting  end  of 
the  resonance  chamber.  The  latter  stagnation  pressure  probe  was  located 
immediately  upstream  of  the  ejector  primary  flow  inlet. 

The  oscillator  knife-edge  can  be  moved  along  its  centerline,  see 
Fig.  3.1-4,  so  that  the  distance  between  the  oscillator  inlet  and  the 
knife-edge  can  be  varied.  However,  the  oscillator  knife-edge  was  always 
positioned  for  the  maximum  oscillator  inlet  to  knife-edge  distance  allowed 
by  the  geometry  of  the  oscillator  body.  Hydraulic  analogy  findings  indi¬ 
cated  that  this  position,  for  the  given  fixed  oscillator  inlet  and  exit 
flow  areas,  would  probably  produce  the  largest  amplitude  pulses  [42], 

This  knife-edge  position  avoided  possible  vibration  problems  which  an 
extended  knife-edge  could  have  created.  However,  this  oscillator  configu¬ 
ration  was  limited  to  only  one  pulse  amplitude  for  a  given  operating 
frequency  and  oscillator  inlet  stagnation  pressure. 

A  planar  two-dimensional,  constant  area,  supersonic-subsonic  ejector 
was  built  for  this  investigation.  The  primary  flow  enters  the  ejector 
mixing  duct  through  a  continuous  slope  converging-diverging  nozzle  designed 
for  a  uniform  Mach  2  exit  plane  flow  by  the  method  of  characteristics  with 
a  boundary  layer  displacement  thickness  correction.  From  a  stagnation 
pressure  traverse  across  the  exit  plane  of  this  nozzle,  a  reasonably  uni¬ 
form  flow  was  found  with  Mach  numbers  in  the  range:  1.98  <  M  <  2.01.  This 
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was  a  pleasant  surprise  because  the  nozzles  were  cut  rather  thin  near  the 
exit  plane  to  minimize  the  base  area  between  primary  and  secondary  flows; 
see  Figs.  3.1-3  and  3.1-6.  The  machining  process  left  the  nozzle  blocks 
slightly  warped  and  it  was  necessary  to  straighten  them.  To  correct 
sealing  problems,  multiple  assembly-disassembly  operations  were  initially 
required.  Based  on  static  pressure  measurements  taken  after  the  final 
assembly  of  the  ejector,  the  centerline  exit  plane  Mach  number  was  found 
to  be  approximately  M  =  1.96. 

The  ejector  unit  itself  consists  of  two  types  of  assemblies.  A  base 
assembly  in  which  secondary  and  primary  flows  are  brought  to  their  con¬ 
fluence  point,  Fig.  3.1-6,  and  three  aluminum  mixing  duct  assemblies, 

Fig.  3.1-7.  The  base  assembly  was  made  primarily  from  stainless  steel  to 
improve  its  dimensional  stability,  to  increase  its  mass,  and  to  reduce 

vibration.  The  supersonic  primary  flow  enters  at  the  center  of  the  mixing 

duct  inlet  interleaved  between  two  subsonic  secondary  streams.  The  nomi¬ 
nal  width,  W,  of  the  constant  area  rectangular  mixing  duct;  i.e.,  the  dis¬ 

tance  across  both  secondary  and  the  primary  streams,  is  33.07  mm.  The 
distance  from  the  primary  flow  centerline  to  a  secondary  flow  wall,  half 
of  the  width,  is  16.54  mm  and  is  indicated  in  Fig.  3.1-6.  The  height  of 
the  mixing  duct,  the  dimension  normal  to  the  page  in  Fig.  3.1-6,  is 
15.24  mm.  The  primary  to  mixing  duct  flow  area  ratio  is  1/3.  Static 
pressure  taps  are  located  along  the  centerline  of  each  of  the  four  mixing 
duct  walls.  The  first  four  pressure  tap  locations  are  at  intervals  of 
X/W  =  1/3  the  duct  width,  W  =  33.07  mm.  The  remainder  of  the  taps  are  at 
one  duct  width  intervals  in  the  X  or  flow  direction  starting  at  X/W  =  2.0. 
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The  base  assembly  pressure  taps  and  the  secondary  stagnation  pressure 
probes  can  be  seen  in  Figs.  3.1-6  and  3.1-3.  Static  pressure  taps  are 
located  in  the  center  of  the  primary  and  secondary  flow  inlets  at  Station 
1  and  on  the  primary  flow  centerline  of  each  mixing  duct  assembly  very 
near  the  exit  plane.  The  three  mixing  duct  assemblies  produce  total  duct 
length- to-width  ratios  of  5,  9,  and  13,  when  mounted  on  the  base  assembly. 
The  L/W  =  5  combination  was  found  to  be  too  short  to  be  of  practical 
interest.  The  shortest  duct  assembly  was  modified  so  that  it  could  be 
attached  to  the  longest  mixing  duct.  Thus,  a  mixing  duct  of  L/W  =  16.374 
could  be  investigated. 

A  back  pressure  valve  and  muffler  were  connected  to  the  exit  plane  of 
the  ejector.  The  back  pressure  valve  was  always  kept  open  and  was  not 
used  to  adjust  the  exit  plane  pressure.  The  valve  was  mounted  on  a  seal¬ 
ing  plate,  thus  providing  a  chamber  into  which  the  ejector  flow  discharged. 
The  muffler  attached  directly  to  the  back  pressure  valve.  A  stagnation 
pressure  probe  was  mounted  on  the  sealing  plate  and  was  used  to  measure 
the  stagnation  pressure  at  the  center  of  the  exit  plane  and  at  a  point  3/8 
of  the  duct  width  off  of  center. 

Figures  3.1-8  and  3.1-9  are  two  views  of  the  assembled  flow  system. 

The  oscillator  is  located  on  the  right  in  both  figures;  the  resonance  and 
feedback  tubes  were  set-up  for  a  142  Hz  operating  frequency.  The  back 
pressure  valve  and  muffler  are  located  on  the  left.  The  holding  tank  is 
behind  the  oscillator. 

An  electronic  data  acquisition  system  capable  of  recording  four 
channels  simultaneously  was  used  to  collect  the  data  for  all  unsteady 
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experiments.  Pressures  were  measured  with  CEC  type  4-312  strain  gage 
pressure  transducers.  These  transducers  were  mounted  in  case  assembles 
or  line  adapters.  Scanivalve  "zero"  volume  line  adapters  were  used; 
these  line  adapters  were  designed  for  use  with  CEC  type  4-312  transducers 
and  are  designed  to  minimize  the  volume  around  the  transducer  diaphragm 
within  the  adapter.  Figure  3.1-10  is  a  cross-sectional  view  of  the 
Scanivalve  "zero"  volume  line  adapter  with  a  transducer  in  place.  Each 
transducer  had  a  cap  custom  fitted  to  its  diaphragm  in  accordance  with 
Scanivalve  specifications.  Pressure  waveform  distortion  due  to  adapter 
and  pressure  line  capacitance  was  minimized  by  using  the  Scanivalve 
adapters  and  by  keeping  the  connecting  pressure  lines  as  short  as  possible. 
A  20  mm  long  piece  of  #16  gage  flexible  tubing  was  used  to  connect  these 
transducer  line  adapters  to  the  pressure  taps  on  the  ejector  unit.  All 
distances  between  ejector  pressure  taps  and  transducer  diaphragms  were 
less  than  75  mm.  CEC  type  1-183  signal  conditioners  were  used  to  excite, 
balance,  and  calibrate  the  transducers  as  well  as  to  provide  an  amplified 
output.  Figure  3.1-11  is  a  schematic  of  the  instrumentation  system  used. 

The  signal  conditioner  output  was  recorded  with  a  Biomation  1015 
waveform  recorder.  This  unit  is  capable  of  recording  four  channels 
simultaneously  at  sampling  rates  as  high  as  100  kHz.  Both  analog  and 
digital  output  were  available  from  the  Biomation,  thus  allowing  the  exami¬ 
nation  of  data  on  an  oscilloscope  before  initiating  the  recording  of  the 
digital  output  on  tape.  The  digital  output  from  the  Biomation  was 
recorded  with  a  Memodyne  model  2146  digital  cassette  recorder.  A  Datos 
305  data  interface  was  used  to  present  the  Biomation  digital  output  in  an 
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acceptable  form  with  formatting  to  the  Memodyne  recorder.  A  Data  Works 
Model  1700  was  used  to  do  this  for  some  of  the  initial  experiments. 

A  Texas  Instruments  Silent  700  ASR  terminal  was  used  to  transmit 
the  recorded  data  to  a  CDC  CYBER  175  computer. 

3.2  EXPERIMENTAL  PROCEDURE 

Cold  flow  air  experiments  were  conducted  to  study  and  compare 
steady  and  unsteady  ejector  performance  characteristics.  All  experiments 
were  conducted  at  stagnation  pressure  ratios,  Ps0/pp0>  that  would  produce 
supersonic  regime  flow  at  sufficiently  low  values  of  P^/P, 0  •  As  can  be 
seen  in  Figs.  2.1-2  and  2.5-4,  the  maximum  pressure  recovery  in  the  back 
pressure  independent  flow  regimes,  which  occurs  at  the  break-off  point 
with  MR  flow,  is  noticeably  less  for  SSR  and  APPR  flows  than  SR  flow. 

Since  this  investigation  is  concerned  with  high  pressure  recovery  applica¬ 
tions  of  ejector  systems,  the  objectives  of  this  investigation  can  most 
efficiently  be  met  by  investigating  only  the  SR  flow  range  of  values  of 
Ps0/Pp0.  Also,  a  pressurized  secondary  system  would  have  been  required  to 
investigate  SSR  and  APPR  flows.  It  can  also  be  seen  in  Fig.  2.1-2  that 
higher  P.^/Pl0»  according  to  the  theoretical  analysis,  are  obtained  in  MR 
flow  than  SR  flow.  Thus  MR  flow,  near  break-off  with  SR  flow,  was  investi¬ 
gated  to  obtain  a  more  complete  understanding  of  ejector  performance  near 
break-off. 

The  effects  of  a  pulsing  driver  flow,  its  frequency,  the  back  pres¬ 
sure,  and  the  mixing  duct  length  on  ejector  performance  and  the  behavior 
of  the  flow  within  the  mixing  duct  were  investigated  for  evidence  of 
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significant  variations  in  mixing  and  pressure  recovery  processes  as  well  as 
ejector  entrainment  characteristics .  Ejector  flow  data  were  also  taken  to 
determine  how  well  the  analysis  predicts  flow  conditions  at  Station  1  for 
both  steady  and  periodic  primary  flow. 

3.2.1  Steady  Flow  Experimental  Procedure 

An  extensive  set  of  steady  flow  experiments  was  conducted  in 
order  to  gain  a  detailed  knowledge  of  the  ejector  operating  characteristics 
and  to  form  a  basis  with  which  an  unsteady  driver  flow  could  be  judged. 

Preliminary  experiments  were  conducted  to  determine  at  what  values  of 
P  ,/P  „  SR  flow  began  to  break-off  to  MR  flow  and  how  w  /w  decreased  as 
Pm3^Pso  was  increased  beyond  the  break-off  value.  These  experiments  were 
accomplished  by  adjusting  the  stagnation  pressures  to  different  levels  so 
tnat  their  ratio,  Ps0/Pp0.  remained  constant  while  the  back  pressure  condi- 
tion,  P^/Pjq.  varied.  During  these  experiments  P^  remained  approximately 
constant.  The  mass  flowrate  ratio  proved  to  be  quite  sensitive  to  the  back 
pressure;  it  decreased  as  pri3/p,  0  increased  once  the  ejector  flow  had 
passed  break-off  conditions  and  was  in  the  mixed  regime.  It  was  learned 
that  the  mixing  duct  static  pressure  profile  could  be  used  as  a  reliable 
indicator  of  whether  or  not  break-off  of  SR  to  MR  flow  had  occurred  and  to 
what  extent.  This  was  helpful  in  conducting  the  experiments  since  the  flow 
began  to  break-off  at  lower  values  of  pm3/ps0  and  more  gradually  than  pre¬ 
dicted  from  the  theoretical  analysis. 

Back  pressure  independent  flow  conditions  were  investigated  for  values 
of  P  /P  in  the  approximate  range:  .048  £  P  /P  <_  .180.  The  mixing 
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duct  static  pressure  distribution  was  used  to  insure  that  the  ejector  flow 
was  in  the  supersonic  flow  regime.  The  pressures:  P  ,  ,  P,,  P  .  ,  P  . , 
and  the  mixing  duct  static  pressure  distribution  were  recorded  by  photo¬ 
graphing  the  manometer  board.  A  precision  pressure  gage  was  used  to  mea¬ 
sure  the  primary  and  exit  plane  stagnation  pressures.  The  necessary  VDI 
nozzle  data  were  also  recorded  to  determine  the  mass  flowrates.  This 
entire  set  of  experiments  was  conducted  for  mixing  ducts  with  length-to- 
width  ratios,  L/W  =  9.0,  13.0,  and  16.374. 

Static  pressure  measurements  versus  time  were  recorded  to  establish 
approximately  the  intensity  and  scale  of  the  turbulence  within  the  mixing 
duct.  As  a  result  of  the  unexpected  nature  of  these  preliminary  transient 
measurements,  a  series  of  experiments  was  conducted  following  the  pro¬ 
cedures  planned  for  the  periodic  driver  flow  experiments.  The  static 
pressures  were  recorded  simultaneously  at  pressure  taps  on  the  primary  flow 
centerline  and  on  the  subsonic  secondary  flow  wall  along  the  length  or  the 
mixing  duct.  These  data  sets  were  recorded  at  values  of  P  /P  „  =  .070  and 
.160  in  the  L/W  =  13.0  mixino  duct  and  at  values  of  P  „/P  n  =  .100  and  .160 
in  the  L/W  =  9.0  mixing  duct.  The  exit  plane  stagnation  pressures  were 
also  recorded.  Static  pressures  within  both  the  L/W  =  13  and  9  mixing 
ducts  at  P  „/P  .  =  .070,  .100,  .130,  and  .160  were  examined  on  an  oscillo- 
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scope  display  which  was  triggered  by  the  Biomation.  The  mixing  duct  flow 
was  examined  at  several  duct  locations  for  side  to  side  symmetry.  The 
secondary  stagnation  pressure,  P  ,  was  examined  for  steadiness. 

Most  of  the  experiments  were  conducted  in  the  mixed  regime  near  the 
experimental  break-off  point  with  the  supersonic  regime.  This  was  done 
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because  the  ejector  flow  began  to  break-off  from  SR  flow  at  values  of 
Pr3/Ps0  that  were  significantly  lower  than  the  theoretical ly  predicted 
values  and  because  the  sharp  theoretical  break-off  between  SR  and  MR  flows 
and  the  rapid  decrease  of  M  in  the  mixed  regime,  as  shown  in  Fig.  2.1-2, 
experimentally  were  found  to  be  a  smoother  transition  followed  by  a  more 
gradual  decline  of  M  ,  and  w  .  Thus,  truly  back  pressure  independent  SR 
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flow  mass  flowrate  ratios,  which  were  obtained  at  the  sacrifice  of  pressure 
recovery,  were  not  significantly  greater  than  MR  flow  values  obtained  at 
higher  P  _/P  MR  flow  near  break-off  was  examined  in  order  to  obtain  a 
more  complete  picture  of  ejector  performance. 

Another  set  of  experiments  was  conducted  for  values  of  Ps0/Pp0  =  .070 
.100,  .130,  and  .160  with  L/W  =  9.0  mixing  duct  to  determine  how  the  mass 
flowrate  ratio  decreased  as  P  ,/P  .  was  increased  above  the  break-off 
value.  The  exit  plane  static  pressure,  P^  ,  was  digitally  recorded  and 
then  averaged  to  determine  an  effective  value  of  P  /Pa 0 .  In  all  other 
respects  the  experimental  procedure  was  the  same  as  that  used  in  the  pre¬ 
liminary  runs  discussed  at  the  beginning  of  this  section;  however,  the 
number  of  experiments  conducted  was  more  extensive. 

3.2.2  Periodic  Flow  Experimental  Procedure 

Periodic  driver  flow  ejector  experiments  were  conducted  at 
driver  frequencies  of  142  Hz  and  250  Hz.  More  experiments  were  conducted 
with  the  L/W  =  9.0  mixing  duct  than  the  two  longer  ducts  because  it  was 
desired  to  determine  whether  or  not  ejector  performance  with  a  shorter 
mixing  duct  could  be  improved  by  a  periodic  driver.  The  mixing  duct  with 
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L/W  =  16.374  duct  was  only  investigated  for  the  back  pressure  independent 
mass  flowrate  ratios  with  the  142  Hz  driver. 

Preliminary  experiments  were  conducted  to  correlate  the  stagnation 
pressure  upstream  of  the  oscillator  inlet  to  the  mean  value,  P  .  ,  of  the 

pO 

stagnation  pressure  in  the  resonance  chamber  of  the  oscillator  slightly 
upstream  of  the  ejector  primary  inlet.  With  this  correlation,  Pp0  could 
accurately  be  set  by  adjusting  the  steady  oscillator  inlet  stagnation 
pressure.  This  correlation  was  found  to  be  constant  within  .5%  over  the 
range  of  pressures  of  interest  and  repeatable  within  calibration  accuracy 
Pp0  was  approximately  90%  of  the  oscillator  inlet  stagnation  pressure  at 
both  driver  frequencies.  It  is  through  this  correlation  procedure  that  a 
representative  or  effective  value  of  Pj0/Pp0  is  obtained  for  comparison 
of  periodic  and  steady  driver  ejector  flow  results. 

The  use  of  an  overbar  to  indicate  the  mean  will  be  dropped  with  the 
mean  value  being  understood.  Proper  reference  to  time  is  made  for 
instantaneous  values. 

The  procedures  outlined  for  the  steady  flow  experiments  using  the 
electronic  data  sampling  equipment  were  followed.  This  includes  static 
pressure  measurements  along  the  length  of  the  ducts  and  mass  flowrate 
break-off  experiments.  These  procedures  were  repeated  completely  with 
the  142  Hz  driver;  for  the  250  Hz  driver  flow,  experiments  were  conducted 
with  only  a  mixing  duct  with  L/W  =  9.0. 

Experiments  were  run  at  both  driver  frequencies  measuring  Pp0,  P  , 
and  1  simultaneously  over  the  range  of  P|0/Pp0-  This  was  done  so  that 
the  ratios  P  ,/P  ,  and  P  ,/P  „  could  be  observed  with  time.  P  ,  could 
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also  be  compared  with  PpQ  to  determine  how  oscillator  pressure  pulses  may 
have  been  altered  by  passage  through  the  C-D  nozzle  at  the  primary  flow 
ini et. 

3.3  EXPERIMENTAL  RESULTS 

Steady  driver  ejector  flow  performance  parameters  have  been  compared, 
in  the  following  discussion,  with  the  one-dimensional  theory  and  with 
periodic  driver  ejector  flow  performance  parameters.  The  comparison  with 
theory  has  been  done  to  identify  trends  in  ejector  performance  over  the 
range  of  flow  conditions  and  mixing  duct  lengths  and  to  test  the  theory. 
The  quasi-steady  analysis  cannot  be  expected  to  predict  periodic  driver 
ejector  flow  conditions  better  than  the  one-dimensional  analysis  predicts 
steady  driver  ejector  flow  conditions.  The  effects  of  a  periodic  driver 
on  ejector  flows  are  judged  by  the  comparisons  made  with  steady  driver 
ejector  flows  and  with  the  predictions  of  the  quasi -steady  analysis. 

3.3.1  Steady  Flow  Experimental  Results 

The  back  pressure  independent  performance  parameters  of  the 

steady  flow  ejector  were  found  to  follow  closely  the  trends  predicted  by 

the  one-dimensional  analysis  with  the  exception  of  the  back  pressure 

ratio,  P irj3 / Ps 0  >  at  the  break-off  point.  In  Fig.  3. 3. 1-1,  it  can  be  seen 

that  the  experimental  and  theoretical  values  of  P|t/P  agree  reasonably 

well.  The  value  of  M  used  in  the  theoretical  analysis  was  determined 

pi 

experimentally  from  P  /P  and  the  isentropic  relations  [43].  The  two 
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longer  mixing  ducts  produced  values  of  Ptl/P  which  were  consistently 
below  the  theoretical  curve;  the  shortest  mixing  duct  produced  values 
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which  were  slightly  above  the  theoretical  curve.  This  suggests  that  1 
and  w  should  be  smaller  for  the  shorter  duct.  Figure  3. 3. 1-2  does  indi- 
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cate  that  lower  w  values  occurred  with  the  shorter  L/W  =9.0  duct  over 
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most  of  the  range  of  Pf Q/P  .  However,  the  mass  flowrate  characteristics 
of  the  shortest  mixing  duct  do  not  follow  the  theoretical  curve  as  well 
as  the  two  longer  ducts,  the  experimental  results  are  below  and  above 
the  theoretical  curve.  It  was  also  observed  that  the  shorter  duct  had  to 
be  operated  at  lower  back  pressure  conditions  before  performance  param¬ 
eters  began  to  level  off  as  supersonic  regime  back  pressure  independence 
was  approached.  The  back  pressure  ratio,  Pn3/Ps0,  for  each  of  the 
experiments  are  indicated  in  Fig.  3. 3. 1-3.  It  did  seem  that  lower  values 
of  P^/P  0  might  have  further  improved  the  shortest  duct  mass  flowrate 
performance,  however,  this  would  have  been  at  the  sacrifice  of  the  result¬ 
ing  ejector  pressure  recovery.  This  matter  was  not  pursued  because  mixed 
flow  static-to-stagnation  pressure  ratios  in  the  exit  plane,  P  ,/P 
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below  the  sonic  flow  value,  .528,  were  obtained. 

During  the  course  of  the  steady  flow  experiments,  the  static  pressure 
profiles  observed  in  all  three  mixing  ducts  changed  dramatically  with 
the  back  pressure  ratio,  P  /Pf Q.  At  back  pressures  lower  than  the  SR 
flow  break-off  value,  the  static  pressure  decreased  smoothly  on  the  sub¬ 
sonic  flow  wall  and  along  the  primary  flow  centerline  as  both  secondary 
and  primary  flows  accelerated  after  entering  the  mixing  duct.  The  flow 
would  remain  at  these  low  pressure  levels  to  within  four  duct  widths  of 
the  exit  plane;  a  rapid  pressure  rise  to  the  exit  plane  condition  would 
take  place  in  the  remainder  of  the  duct.  As  P  ,/P  „  was  increased,  the 
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distance  downstream  in  the  duct  that  the  flow  would  remain  in  a  low  pres¬ 
sure,  high  velocity  state  decreased.  The  value  of  Pri5/Pl0  at  which  the 
initial  pressure  decrease  and  the  corresponding  acceleration  was 
immediately  followed  by  the  onset  of  the  pressure  recovery  process  was 
approximately  the  break-off  value  of  P^/P^-  A  decrease  in  the  mass  flow- 
rate  ratio  was  always  observed  to  occur  as  P  ,/P  „  was  further  increased 
for  a  given  value  of  Pj0/Pp0  in  the  L/W  =  13.0  and  L/W  =  16.374  mixing 
ducts.  Thus,  break-off  from  SR  to  MR  flow  had  occurred.  As  P  ,/P  .  was 
further  increased  the  shock  structure  in  the  supersonic  core  became 
stronger,  particularly  at  lower  values  of  Ps0/Pp0-  Although  the  primary 
flow  expanded  and  accelerated  to  low  pressures  as  it  did  at  lower  values 
of  Pm3/Pj0,  the  secondary  flow  was  accelerated  less  as  indicated  by  the 
smaller  pressure  drop  along  the  wall  downstream  of  the  entrance  of  the 
duct.  Station  1.  As  P^j/P, 0  was  increased  further,  secondary  flow  wall 
pressures  in  the  acceleration  section  of  the  duct  remained  constant  and 
equal  to  P  . .  At  this  value  of  P  _/P  n,  the  mass  flowrate  ratio  had 
probably  decreased  by  approximately  50%  from  the  SR  flow  value. 

The  important  features  of  the  mixing  duct  static  pressure  profile  are 
demonstrated  in  Fig.  3. 3. 1-4.  Pressure  taps  along  the  subsonic  secondary 
flow  walls  are  referred  to  as  being  at  the  wall  and  pressure  taps  on  the 
primary  flow  centerline  are  referred  to  as  being  at  the  center  in  all  of 
the  following  discussions  and  in  all  figures.  The  wall  and  center  tap 
static  pressures  along  the  length  of  the  L/W  =  16.374  mixing  duct  are  pre¬ 
sented  in  Fig.  3.3. 1-4  for  two  steady  flows  with,  as  near  as  possible,  the 
same  value  of  Ps0/Pp0-  Although  the  difference  between  the  two  profiles 
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is  significant,  the  profile  for  the  lower  value  of  P^/P^  was  obtained  by 
increasing  the  value  of  PpQ  above  that  of  the  other  profile  from  71.22  psia 
to  73.80  psia  and  adjusting  Ps0  accordingly.  With  a  periodic  driver,  the 
reverse  case  v/as  found,  i.e.,  the  value  of  P^j/P  0  is  approximately  con¬ 
stant  while  Ps0/Pp0  changes;  however,  the  effects  on  the  mixing  duct  pres¬ 
sure  profile  are  still  the  same.  At  these  pressure  levels,  the  oscillator 
pulses  have  a  total  amplitude  of  12  to  15.25  psi.  Using  the  theoretical 
results  as  a  guide,  the  effective  changes  in  the  break-off  value  of 
P  ,/P  „  should  be  larger  than  the  difference  between  the  two  values  shown 
in  Fig.  3. 3. 1-4  at  the  sjme  value  of  Ps0/Pp0- 

3.3.2  Periodic  Flow  Experimental  Results 

During  the  experiments  with  the  steady  and  periodic  driver 
flows,  three  types  of  unsteady  flow  phenomena  were  observed  in  the  mixing 
duct  flow.  The  largest  disturbances  found  within  the  mixing  duct  appear 
to  be  due  to  changes  in  and  motion  of  the  shock  structure  in  the  super¬ 
sonic  core.  These  shock  waves  and  their  location  seemed  quite  unstable, 
which  may  have  been  due  to  other  unsteady  disturbances.  A  low  amplitude 
oscillation  of  the  static  pressures  within  the  duct  with  frequencies  up  to 
approximately  1000  Hz  was  also  observed.  This  oscillation  was  more 
noticeable  in  the  mixing  duct  with  L/W  =  9.0  than  with  L/W  =  13.0.  These 
oscillations  were  present  in  the  wall  or  secondary  flow  at  Station  1  but 
if  observed,  were  not  seen  in  the  center  flow  until  after  X/W  =  1.0. 
Oscillation  amplitudes  usually  increased  in  the  acceleration  section  of 
the  mixing  duct;  however,  low  back  pressure  conditions  were  found  to 
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inhibit  this  growth  and  to  diminish  the  amplitudes  at  Station  1.  The  most 
surprising  phenomena  found  with  either  steady  or  periodic  driver  flows  are 
shown  in  Fig.  3. 3. 2-1.  Large  amplitude  pulses,  sometimes  bursting  into 
periodic  groups,  were  observed  in  most  flows.  These  pulses  originated  in 
the  wall  flow  at  approximately  the  end  of  the  acceleration  section  of  the 
flow  which  is  near  the  location  of  the  first  shocks  embedded  in  the  center 
flow.  The  amplitude  and  frequency  of  occurrence  of  these  bursts  increased 
as  Pj0/Pp0  decreased.  In  many  cases  these  bursts  appear  to  have  been  an 
amplification  of  the  smaller  amplitude  pulses  mentioned  above.  Neither 
the  origin  or  amplification  of  these  disturbances  is  well  understood. 
Changes  in  operating  conditions  which  adversely  increased  the  secondary 
flow  pressure  gradient  through  the  duct  and  increased  the  strength  of 
the  shocks  in  the  supersonic  core  are  believed  to  stimulate  these  pulses. 

It  appears  that  corner  flow  separation  and  reattachments,  shear  layer 
instabilities,  and  disturbances  due  to  the  unstable  shock  structure  which 
propagate  into  the  subsonic  wall  flow  are  probable  mechanisms  for  these 
unsteady  phenomena.  Investigations  of  unsteady  turbulent  boundary  layers 
have  found  bursts  in  intensity  levels  do  occur  in  some  cases  with  periodic 
freestream  flows  [44,45].  The  turbulence  at  the  edge  of  the  boundary 
layer  can  "violently"  increase  in  level  and  burst  into  the  freestream  flow. 
Oscillations  in  the  shock  structure  of  the  supersonic  core  may  be  coupled 
with  larger  scale  turbulence  to  stimulate  such  bursts.  Figure  3. 3. 2-2 
shows  one  more  feature  of  all  pulses,  that  is  side-to-side  symmetry  was 
always  found  in  all  flows. 
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A  set  of  static  pressure  profile  data  taken  along  the  length  of  the 
mixing  duct  with  L/W  =9.0  and  with  Pj0/Pp0  =  -100  for  steady  and  both 
periodic  driver  flows  is  included  in  Figs.  3. 3. 2-3  to  3.3.2-21  to  present 
a  clear  picture  of  what  occurs  within  the  mixing  duct  and  how  some  of  the 
many  variables  affect  the  flow.  All  data  was  taken  for  steady,  142  Hz, 
and  250  Hz  driver  flows  with  approximately  equivalent  input  and  back  pres¬ 
sure  flow  conditions.  Figures  3. 3. 2-3  and  3. 3. 2-4  are  the  primary  and 
secondary  static  pressures  at  Station  1  versus  time  for  both  periodic 
driver  frequencies.  The  secondary  flow  pulses  are  somewhat  irregular  and 
the  periodicity  appears  to  be  much  weaker  at  the  lower  frequency.  These 
irregular  oscillations  in  the  secondary  flow  were  not  as  large  in  the  mix¬ 
ing  duct  with  L/W  =  13.0.  The  oscillations  in  the  secondary  flow  have 
grown  by  X/W  =  0.333  and  are  prominent  with  both  steady  and  periodic 
driver  flows;  these  data  are  shown  in  Figs.  3. 3. 2-5,  3. 3. 2-6,  and  3. 3. 2-7. 
Note  that  the  center  flow  has  not  been  affected  by  the  irregular  disturb¬ 
ances  observed  in  the  subsonic  wall  flow.  Also,  note  that  the  250  Hz  wall 
flow  pulses  slightly  precede  those  in  flow  at  the  center.  This  phase 
shift  can  be  seen  in  nearly  all  of  the  data  taken  in  the  first  1/3  of  the 
mixing  duct;  it  was  usually  more  noticeable  with  the  250  Hz  driver  flow. 
Static  pressures  continue  to  drop  as  the  flows  interact  and  accelerate 
until  somewhere  near  X/W  =  1.0.  Shocks  occur  in  the  primary  flow  where  it 
stops  expanding  against  the  secondary  and  readjusts  its  flow  direction. 
Disturbances  that  occur  in  the  wall  flow  also  appear  in  the  center  flow  at 
X/W  =  1.0,  see  Figs.  3. 3. 2-8,  3. 3. 2-9,  and  3.3.2-10.  The  periodicity  of 
the  center  flow  seems  to  fade  in  and  out.  The  pulses  at  X/W  =  1.0  in  the 
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expanded  142  Hz  driver  ejector  flow  are  larger  and  steeper  than  could  be 
produced  solely  by  driver  flow  pulsations.  Figure  3.3.2-11  is  an  example 
of  how  much  the  driver  pulses  were  diminished  by  the  expansion  of  the 
supersonic  primary  flow.  The  data  in  this  figure  were  taken  for  a  back 
pressure  independent  SR  flow  experiment  with  Ps0/Pp0  -  .070.  The 
Ps0/Pp0  =  -100  data  under  consideration  has  a  back  pressure  sufficiently 
high  for  the  ma;,s  flowrate  ratio  to  have  decreased  about  9%  from  the  back 
pressure  independent  SR  flow  value.  Shock  waves,  due  to  the  relatively 
high  back  pressure,  are  present  in  this  MR  flow  and  their  unsteady  behavior 
in  the  periodic  flow  is  believed  to  be  the  cause  of  the  larger  pulses  in 
the  flow  at  X/W  =  1.0.  As  mentioned  earlier,  the  strength  and  in  most 
cases  the  number  of  shocks  in  the  flow  increase  when  SR  flow  breaks-off 
to  MR  flow  due  to  an  increase  in  P  ,/P  „ .  Note  also  that  the  unsteadiness 
of  the  shock  waves  in  the  steady  driver  ejector  flow  is  comparable  to  the 
periodic  driver  ejector  flows.  Shock  waves  are  clearly  present  at  the 
X/W  =  1.333  pressure  tap  location,  as  shown  in  Figs.  3.3.2-12,  3.3.2-13, 
and  3.3.2-14.  The  smaller  oscillations  at  the  subsonic  wall  flow  pressure 
taps  are  similar  in  appearance  to  the  large  oscillations  at  the  center 
taps  and  are  phase-shifted  slightly  ahead  of  the  disturbances  at  the 
center.  The  periodic  driver  moves  the  shock  waves  upstream  and  downstream 
within  the  mixing  duct;  some  of  these  shock  waves  may  disappear  or 
coalesce  and  then  reappear.  The  shock  waves  in  the  steady  driver  ejector 
flow  are  quite  unsteady;  however,  their  unsteady  action  is  intermittent 
and  not  as  vigorous  as  in  the  periodic  driver  flows.  The  shock  wave 
patterns  are  the  dominant  feature  within  the  flow  in  the  mixing  duct  until 


X/W  =  6.0;  however,  at  X/W  =  4.0,  the  wall  flow  appears  to  be  decoupled  from 
the  action  in  the  primary  flow  at  the  center.  These  mixing  duct  flow 
features  are  shown  in  Figs.  3.3.2-15,  3.3.2-16,  and  3.3.2-17.  The  250  Hz 
driver  flow  again  appears  more  pulsatile  than  the  142  Hz  driver  flow.  The 
static  pressure  levels  at  X/W  =  4.0  have  not  increased  significantly  since 
X/W  =  1.333  and  the  flow  does  not  appear  to  have  mixed  significantly 
either.  Shock  waves  are  present  in  the  flow  downstream  of  X/W  =  6.0  but 
are  not  evident  at  X/W  =  8.0.  At  X/W  =  8.0,  the  center  and  wall  flows  have 
become  noticeably  more  broken  and  turbulent.  The  mixing  duct  flow  is  sig¬ 
nificantly  more  uniform  across  the  duct;  the  mixing  duct  flow  at  this  loca¬ 
tion  is  essentially  the  same  for  all  three  of  the  driver  flows,  see 
Figs.  3.3.2-13,  3.3.2-19,  and  3.3.2-20.  The  142  Hz  driver  ejector  flow  is 
weakly  periodic  at  X/W  =  8.0  while  the  250  Hz  driver  ejector  flow  is 
definitely  periodic.  The  supersonic  core  flow  seems  to  have  finally  disin¬ 
tegrated  in  all  three  driver  flows  between  X/W  =  6.0  and  X/W  =  8.0,  thus 
initiating  a  rapid  mixing  and  pressure  recovery  process.  This  rapid  mixing, 
which  occurred  immediately  downstream  of  the  disappearance  of  the  supersonic 
core,  was  also  noticed  in  water  table  studies  of  supersonic-subsonic  ejector 
flow  where  dye  was  injected  into  the  wall  flow  to  indicate  mixing  rates. 

At  the  exit  plane,  all  three  driver  ejector  flows  are  nearly  identical. 
Periodicity  was  still  weakly  present  at  the  exit  plane  of  the  mixing  duct 
due  to  the  muffler  and  back  pressure  valve  flow  resistance;  this  is  shown 
in  Fig.  3.3.2-21.  It  should  be  noted  that  the  mean  static  pressure 
levels  at  most  locations  in  the  mixing  duct  are  approximately  the  same  for 
all  three  driver  flows.  As  shown  in  Figs.  3.3.2-22  and  3.3.2-23,  the  exit 


56 


plane  stagnation  pressures  were  nearly  indistinguishable  and  periodicity 
was  not  a  strong  feature. 

The  only  notable  difference  between  the  steady  and  periodic  driver 
ejector  flows  in  the  mixed  regime  flows  discussed  above  was  the  increased 
activity  of  the  shock  waves  in  the  supersonic  primary  flow.  Supersonic 
regime  flows  operating  at  the  same  Ps0/pp0  as  a  mixed  regime  flow  had 
fewer  and  weaker  shocks  than  the  mixed  regime  flow.  Thus, supersonic  regime 
flows  were  less  susceptible  to  disturbances  due  to  a  periodic  driver  flow 
than  mixed  regime  flows. 

The  primary-to-secondary  stagnation  pressure  ratio  had  a  pronounced 
effect  on  mixing  duct  flow  features.  Lower  values  of  the  stagnation  pres¬ 
sure  ratio  tended  to  have  strong  shock  waves  within  the  flow  for  all  but 
the  low  back  pressure  SR  flows.  The  experimental  data  with  Ps0/Pp0  =  .100 
that  was  just  discussed  is  typical. 

The  opposite  case  is  true  for  higher  values  of  ps0/pp0-  For 
P  _/P  .  =  .160,  the  resultant  P  ,/P  ,  is  roughly  .9.  In  this  case  the 

s  0  p  0  si*  pi  Zf  u 

primary  flow  does  not  expand  as  much  after  entering  the  mixing  duct  as  it 
did  at  the  lower  values  of  Ps0/Pp0;  consequently,  the  resulting  shock 
waves  in  the  supersonic  core  flow  do  not  have  as  great  an  effect  on  the 
flow.  Static  pressure  measurements  were  taken  in  the  L/W  =  9.0  mixing  duct 
with  Ps0/Pp0  =  -160  and  a  back  pressure  ratio  sufficiently  high  to  have 
caused  the  mass  flowrate  ratio  to  decrease  by  approximately  8%  from  the 
back  pressure  independent  value.  The  largest  disturbances  within  this 
flow  are  shown  in  Figs.  3.3.2-24,  3.3.2-25,  and  3. 3. 2-^6.  The  pulses  in 
the  periodic  flow  were  probably  due  to  a  weak  shock  wave  moving  and/or 
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changing  in  strength.  This  was  the  only  location  where  disturbances  of 
any  notable  size  were  found.  The  static  pressures  at  Station  1  for 
P  /P  =  .160  and  the  L/W  =  9.0  mixing  duct  and  for  P  /P  =  .070  and 

sOpO  sUpU 

the  L/W  =  13.0  mixing  duct  are  shown  in  Figs.  3.3.2-27  and  3.3.2-28 
respectively.  From  these  figures,  the  secondary  flow  pulses  are  seen  to 
decrease  appreciably  with  Ps0/Pp0- 

Experimental  values  of  Psl/Ppl  versus  time  were  compared  with  values 
determined  from  the  theoretical  quasi-steady  analysis.  Obtaining  the 
P  /P  „  versus  time  data  as  required  for  the  theoretical  analysis  was  not 

s  0  pO 

as  straightforward  as  was  initially  expected.  Figures  3.3.2-29  and 
3.3.2-30  are  the  simultaneously  recorded  primary  stagnation  and  static 
pressures  for  the  142  Hz  and  250  Hz  drivers,  respectively.  The  stagnation 
pressure  probe  was  located  immediately  upstream  of  the  primary  C-D  nozzle 
inlet  and  the  static  tap  was  located  at  the  center  of  the  nozzle  exit 
plane.  Station  1.  The  separation  distance  between  these  two  locations  led 
to  the  small  phase  shift  seen  in  both  figures.  The  250  Hz  driver  flow 
stagnation  pressure  profile  was  pressure  invariant.  The  higher  frequency 
components  in  the  142  Hz  driver  flow  profiles  grew  slightly  in  amplitude, 
relative  to  that  of  the  base  frequency,  as  the  mean  pressure  level 
increased.  Difficulties  were  encountered  when  examining  the  variation 
of  P  ,/P  „  with  time.  The  decision  was  made  to  phase-shift  P  ,  with 
respect  to  P  „  such  that  the  relatively  sharp  valleys  of  each  cycle  were 

p  0 

in  phase  before  determining  P  /P  as  a  function  of  time.  P  /P  was 

p 1  p0  p 1  pO 

found  to  be  oscillatory  and  to  be  affected  strongly  by  the  higher  frequency 
components  of  P  which  were  not  present  in  P  ;  this  is  shown  in 
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Figs.  3.3.2-31  and  3.3.2-32.  Note  that  for  the  250  Hz  driver  flow,  the 
P  ,  /P  .  curve  was  determined  from  a  different  set  of  data  than  was  used  in 

pi  pO 

Fig.  3.3.2-30;  both  sets  of  data  were  in  phase  with  one  another  and  the 
pressure  difference  between  the  two  was  not  large.  The  same  data  was  used 
to  determine  both  142  Hz  plots.  The  ratio,  Ppl/Pp0»  for  the  250  Hz  driver 
flow  was  at  a  minimum  when  P  „  was  at  a  maximum  but  the  142  Hz  data  seems 

pO 

to  indicate  just  the  opposite.  The  mean  values  of  P  ,/P  „  at  both  driver 

p  1  pO 

frequencies  was  approximately  constant  at  a  value  of  P  , /P  „  =  .140  as  com- 

p  1  pO 

pared  to  the  steady  flow  value  of  Ppl/Pp0  =  .136.  The  pulsed  flow  in  the 
resonance  chamber  of  the  oscillator  adjusts,  by  wave  action,  to  the  area 
change  between  the  resonance  chamber  and  C-D  nozzle  throat  as  well  as  the 
choked  flow  condition  at  the  nozzle  throat  [46].  This  wave  action  adjust¬ 
ment  of  the  flow  apparently  has  altered  the  wave  form  of  the  pressure 
pulses  between  the  stagnation  probe  and  Station  1  as  indicated  in 
Figs.  3.3.2-29  and  3.3.2-30.  In  order  to  obtain  P  _/P  n  versus  time  for 
the  theoretical  analysis,  values  of  P  were  recorded  and  divided  by  the 
mean  value  of  P  ,/P  „  at  the  given  pressure  conditions  to  determine  P  „ 

pi  p  O  ar  p0 

versus  time.  P  „  was  held  constant  and  P  ,  P  ,  ,  and  P  „  were  all  simulta- 

s  O  s  1  p 1  ’  pO 

neously  measured.  The  agreement  between  experimental  and  theoretical 
results  is  reasonably  good;  these  comparisons  are  in  Fig.  3.3.2-33.  The 
experimental  data  with  a  mean  value  of  P  ,/P  ,  =  .504  was  taken  for  a  flow 

S  1  pi 

in  the  mixed  regime  with  the  back  pressure  sufficiently  high  to  have  caused 
an  8%  decrease  in  the  mass  flowrate  ratio  from  the  supersonic  regime  value. 
Thus,  M  and  w  /w  are  less  than  the  break-off  value  while  P  /P  is 

s  1  s  p  slpl 

greater.  The  data  with  a  mean  value  of  P  /P  =  .656  was  taken  for  a  flow 
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in  the  supersonic  regime  while  the  data  with  a  mean  value  of  Pfl/P  =  .875 
was  taken  for  a  flow  in  the  mixed  regime  with  a  mass  flowrate  3  to  4%  less 
than  the  break-off  value.  The  value  of  the  primary  nozzle  exit  Mach  number 
used  in  the  analysis  was  M  =  1.957,  the  same  value  as  was  used  in  the 
steady  flow  analysis. 

The  imposed  back  pressure  had  other  effects  on  the  performance  of  the 
ejector  than  those  previously  described.  Depending  on  the  back  pressure, 
the  ejector  exit  plane  stagnation  pressures  were  found  to  differ  between 
the  wall  and  center  probe  locations  by  as  much  as  a  factor  of  two  or  by  very 
little.  Lower  back  pressures  and  the  shortest  mixing  duct  produced  the 
most  non-uniform  flow  at  the  ejector  exit  plane.  The  periodic  driver  flow 
had  little  effect  on  the  exit  plane  stagnation  pressure  uniformity.  The 
mean  value  of  the  center  and  wall  stagnation  pressures  were  more  nearly 
equal  for  some  periodic  flows  that  were  near  the  break-off  back  pressure 
ratio.  However,  results  were  inconsistent.  The  ratio,  Ps0/Pp0>  seemed  to 
have  an  effect  and  Pni3/Ps0  itself  produced  a  uniform  exit  plane  flow  if  it 
was  sufficiently  high. 

As  was  mentioned  earlier,  SR  ejector  flows  at  low  back  pressures 
remained  in  a  low  pressure  high  velocity  state,  free  of  shock  waves  of 
noticeable  strength,  to  within  four  duct  widths  of  the  exit  plane.  It  was 
observed  that  under  these  flow  conditions  the  wall  and  center  pressure 
pulses  were  shifted  completely  out  of  phase  downstream  of  the  initial 
acceleration  section  of  the  mixing  duct;  this  is  shown  in  Fig.  3.3.2-34. 

An  expansion  of  the  supersonic  primary  flow,  as  indicated  by  a  pressure 
drop  at  the  center  pressure  tap,  caused  a  deceleration  and  pressure  rise 
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in  the  secondary  flow  at  the  wall  pressure  tap.  Thus,  it  appears  that 
after  choking  the  secondary  flow  accelerates  to  a  supersonic  flow  if  the 
primary  contracts  slightly  due  to  shock  waves  imbedded  in  the  supersonic 
core.  Nelson  [47]  found  that  normal  shock  waves  occurred  in  the  wall  flow 
of  a  similar  two-dimensional  constant  area  supersonic-subsonic  ejector. 

The  way  that  the  mass  flowrate  ratio  declines  when  an  SR  flow  breaks- 
off  to  MR  flow  due  to  an  increase  in  P  ,/P  „  was  found  to  differ  for 
steady,  142  Hz,  and  250  Hz  driver  flows  at  most  values  of  P,0/Pp0  which 
were  examined.  It  was  observed  that  the  back  pressure  independent  SR  mass 
flowrate  ratios  were  approximately  the  same  for  all  three  driver  flows  but 
that  the  MR  flow  mass  flowrate  ratios  at  a  given  back  pressure  ratio, 
Pmj/Pjo’  were  higher  with  the  periodic  drivers  than  with  the  steady  driver 
Figures  3.3.2-35,  3.3.2-36,  and  3.3.2-37  show  these  mass  flowrate  trends. 
It  appears  that  the  periodic  driver  ejector  flows  break-off  at  higher 
P  /Pj0  than  the  steady  driver  ejector  flows.  The  maximum  difference 
between  the  values  of  P  ,/P  „  at  a  given  P  /P  .  and  the  w  /w  occurs 
between  the  steady  and  250  Hz  driver  ejector  flows  at  Ps0/Pp0  =  .070;  this 
difference  is  less  than  5%.  The  142  Hz  driver  did  not  differ  from  the 
steady  driver  by  a  significant  amount  in  most  cases  and  all  three  driver 
ejector  flows  followed  the  same  break-off  trends  at  P  /P  =  .160. 
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Figure  3.1-1  Flow  schematic  of  the  test  rig 
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Figure  3.1-11  Schematic  detailing  the  data  acquisition  system 
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pressure  independent  steady  flow  values  of  P  /P  versus 


pressure  independent  steady  flow  values  of  w  /w  versus 
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Figure  3. 3. 2-2  Wall  flow  static  pressures  were  found  to  be  symmetric  about  the 
primary  flow  centerline  (L/W  =  9.0,  P  /P  -  .100) 


Figure  3. 3. 2-4  Primary  and  secondary  static  pressures  at  Station  1  for  the  250  Hz  driver 
ejector  flow  versus  time  (L/W  =  9.0,  P  /P  =  .100) 
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pressures  versus  time  in  the  steady  driver  ejector  flow 
=  1 .00  (L/W  =  9.0,  P  n/P  =  .100) 
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Figure  3.3.2-11  Static  pressure  pulses  in  a  highly  expanded  shoe 
is  relatively  low  (L/W  =  13.0,  P  /P  =  .070) 
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pressures  versus  time  in  the  steady  driver  ejector 
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pressures  versus  time  in  the  142  Hz  driver  ejector 
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pressures  versus  time  in  the  steady  driver  ejector 


pressures  versu:  time  in  the  250  Hz  driver  ejector 
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Figure  3.3.2-22  Exit  plane  stagnation  pressures  versus 


PLANE  TOTAL  PRESSURES 
vs.  TIME 


Exit  plane  stagnation  pressures  versus  time  in  the  250  Hz  driver 
ejector  flow  (L/W  =  9.0,  P  / P  -  .100) 


pressures  versus  time  in  the  steady  driver  ejector 


PRESSURES  A 
vs.  TIME 


pressures  versus  time  in  the  142  Hz  driver  ejector 


pressures  versus  time  in  the  250  Hz  driver  ejector  flow 


STATIC  PRESSURES  AT  STATION  1 


3.3.2-27  Static  pressures  at  Station  1  versus  time 


Figure  3.3.2-28  Static  pressures  at  Station  1  versus  time 
(L/W  =  13.0,  P  /P  n  =  .070) 


Simultaneously  recorded  values  of  P  t  and  P  versus  time 
at  the  250  Hz  driver  frequency  ’ 


versus  time  at  the  142  Hz  driver  frequency 


Figure  3.3.2-32  P  /P  versus  time  at  the  250  Hz  driver  frequency 


PRESSURE  RATIO  AT  STATION  1 
vs.  TIME 


Theoretical  and  experimental  values  of  P  /P  versus 


pressures  versus  time  at  X/W  =  1.333  with  a  relatively  low 
escure  (L/W  =  13.0,  P  /P  -  .070) 


MASS  FLOWRATE  RATIO  vs.  BACK  PRESSURE 
PBD/Pp(.  =  .070  L/W  =  9.0 


for  flows  that  have  passed  the  break-off 


versus  P  /P  for  flows  that  have  passed  the  break-off 


4.0  CONCLUSIONS 


The  conclusions  that  can  be  drawn  as  a  result  of  this  experimental 
and  theoretical  investigation  of  multiple  ducted  streams  with  a  steady  or 
a  periodic  supersonic  driver  flow  are: 

1.  The  one-dimensional  analysis  predicts  the  back  pressure 
independent  values  of  ws  /wp  and  Psl/P  for 

the  steady  driver  ejector  flow  reasonably  well  over  the 
range  of  values  of  Ps0/Pp0  and  mixing  duct  lengths 
investigated. 

2.  The  values  of  the  back  pressure  ratio,  PI1-3/Ps0>  at  which 
both  steady  and  periodic  driver  ejector  flows  began  to 
break-off  from  back  pressure  independent  SR  flow  to  MR 
flow  were  25  to  35  percent  lower  than  theoretically 
predicted. 

3.  Periodic  primary  ejector  flows  were  found  to  drive  and 
determine  the  secondary  flow  in  such  a  way  that  the 
quasi-steady  assumption  appears  to  be  valid.  Theoretical 
and  experimental  values  of  P^/P  versus  time  were  in 
reasonably  good  agreement.  Back  pressure  independent 
values  of  the  mass  flowrate  ratios  for  the  periodic  driver 
ejector  flows  were  in  good  agreement  with  the  mean  values 
predicted  by  the  quasi-steady  analysis.  The  steady  flow 
values  of  the  mass  flowrate  ratios  at  the  corresponding 
mean  values  of  P  /P  were  nearly  the  same  as  the 

s  0  p  0 

periodic  driver  flow  values. 
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4.  Mass  flowrate  ratios  and  pressure  recovery  capabilities 
of  the  ejector  were  nearly  the  same  for  the  L/W  =  16.374 
and  L/W  =  13.0  mixing  ducts.  The  L/W  =  9.0  mixing  duct 
flow  broke-off  from  SR  to  MR  flow  at  lower  values  of  the 
back  pressure  ratio  than  flow  in  the  longer  two  mixing 
ducts.  Mass  flowrate  ratios  were  less  in  the  shortest 
mixing  duct  at  given  values  of  P  / P  „  and  P  ,/P  „  than 
in  the  two  longer  mixing  ducts. 

5.  Periodic  driver  ejector  flows  subject  to  back  pressures 
high  enough  to  have  caused  break-off  from  SR  to  MR  flow 
operated  at  higher  P  ,/P  ,  up  to  5  percent  higher  than 
the  steady  flow  ejector,  at  given  values  of  P  /P  „  and 

s  0  p  0 

w  /w  in  the  L/W  =  9.0  mixing  duct.  Explained  another 
*  p 

way,  one  might  say  that  higher  mass  flowrate  ratios  were 
obtained  at  given  values  of  P  A/P  .  and  P  ,/P  for 
ejector  flows  that  had  passed  the  break-off  point  and 
were  in  the  mixed  regime.  However,  this  is  not  felt  to 
be  the  result  of  increased  entrainment  due  to  the 
periodicity  of  the  driver  but  rather  due  to  a  shift  in 
the  position,  towards  the  SR  break-off  point,  that  the 
ejector  operates  at  on  the  MR  solution  surface,  see 
Fig.  2.1-2,  for  a  given  value  of  Pj0/Pp0.  The  agreement 
between  steady  and  periodic  driver  flow  back  pressure 
independent  mass  flowrate  ratios,  the  quasi -steady  nature 
of  the  periodic  ejector  flows  at  Station  1,  and  the 
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linearity  of  the  theoretical  values  of  w  /w  as  a  func- 

»  p 

tion  of  Ps0/Pp0>  see  Fig.  3. 3. 1-2,  suggests  that  the 
entrainment  of  the  secondary  flow  in  the  mixed  regime 
near  break-off  with  the  supersonic  regime  is  not  sig¬ 
nificantly  affected  by  the  periodicity  of  the  driver 
flow. 

6.  Static  pressure  measurements  versus  time  did  not  indicate 
that  significant  differences  in  the  pressure  recovery 
process  or  flow  characteristics  within  the  mixing  duct 
existed  between  periodic  and  steady  driver  ejector  flows. 
Shock  waves  in  the  supersonic  primary  flow  were  greatly 
disturbed  by  a  periodic  driver  but  the  shock  waves  in 
the  steady  driver  ejector  flow  were  also  quite  transient 
and  unstable.  Mean  values  of  static  pressures  at  the 
same  location  and  under  corresponding  flow  conditions 
were  approximately  the  same  in  periodic  and  steady  driver 
ejector  flows. 

7.  Mixing  rates  between  the  primary  and  secondary  flows 
were  not  noticeably  altered  by  a  periodic  driver.  The 
location  of  the  point  at  which  the  static  pressures 
became  nearly  uniform  across  the  mixing  duct  was  approx¬ 
imately  the  same  for  periodic  and  steady  driver  ejector 
flows.  Exit  plane  stagnation  pressure  uniformity  was 
slightly  improved  in  only  a  few  cases  by  a  periodic 
driver  flow.  It  appears  that  rapid  mixing  did  not  take 
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place  until  the  primary  flow  in  the  center  of  mixing 
duct  became  subsonic.  As  indicated  by  the  presence  of 
shock  waves,  the  periodicity  of  the  primary  flow  did 
not  change  the  average  location  of  the  disintegration 
of  the  supersonic  primary  to  subsonic  flow.  The 
periodic  driver  had  little  effect  on  the  character  of 
the  flow  in  the  high  mixing  rate  region  which  lies 
between  the  approximate  location  of  the  loss  of  the 
supersonic  primary  flow  and  the  exit  plane  of  the  mixing 
duct. 

8.  The  Strounal  number  of  the  periodic  driver  ejector  flows 
were  low  in  comparison  to  those  in  the  mixing  studies 
cited  in  Section  1.2.  If  based  on  the  velocity  of  the 
primary  flow  at  Station  1,  the  hydraulic  diameter  of  the 
primary  flow  C-D  nozzle,  and  the  250  Hz  driver  frequency, 
the  Strouhal  number  is  St  =  .0064.  The  Strouhal  number 
for  a  subsonic  secondary  stream  at  Station  1  with  the 
theoretical  flow  velocity  at  P  „/P  „  =  .100  and  with 

S  O  p  U 

w  =  250  Hz  is  St  =  .025.  With  the  theoretical  value  of 
Mm3  ^or  Ps0/pPo  =  ‘I 00  and  the  mixing  duct  hydraulic 
diameter,  the  Strouhal  number  at  the  exit  plane  is 
St  =  .029  with  co  =  250  Hz.  A  fluid  particle  will  pass 
entirely  through  the  mixing  duct  in  only  a  fraction  of  a 
driver  cycle.  The  driver  flow  frequencies  used  in  this 
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investigation  were  too  low  to  significantly  effect  the 
dynamics  of  mixing  between  high  speed  compressible 
streams  in  a  snort  mixing  duct. 

9.  The  large-scale  fluidic  oscillator  was  an  efficient  and 
reliable  method  of  producing  periodic  pulses.  It  is  not 
felt  that  frequencies  higher  than  approximately  500  Hz 
could  be  obtained  with  a  fluidic  oscillator  configuration 
similar  to  one  used  in  this  investigation  without  decreas¬ 
ing  the  size  of  the  oscillator  to  such  an  extent  that  the 
mass  flowrate  through  the  device  must  be  appreciably 
decreased. 

10.  The  transient  disturbances  in  the  steady  driver  ejector 
flows  were  comparable  to  disturbances  produced  by  the 
periodic  driver  in  most  cases;  thus,  the  effects  of  a 
periodic  driver  on  ejector  flow  have  not  been  clearly 
delineated  by  the  "steady"  flow  comparison.  An  axisym- 
metric  ejector  configuration  may  prove  to  be  more  stable 
than  2-D  planar  configuration  examined  and  may  therefore 
provide  a  more  meaningful  comparison. 
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APPENDIX  A 

UCAE,  Unsteady  Constant  Area  Ejector  program,  was  written  to  perform 
the  quasi-steady,  one-dimensional  analysis  discussed  in  Section  2.0.  The 
program  and  a  sample  input  and  output  are  listed  in  this  appendix. 

Default  values  of  UCAE  assume  air-to-air  ejector  flow  with  equiva¬ 
lent  primary  and  secondary  stagnation  temperatures.  Variable  definitions 
are  given  at  the  beginning  of  the  program.  MP1 ,  AP1M3,  PM3S0I,  and  MMAX 
are  not  assigned  default  values  and  must  be  assigned  values  in  the  first 
NAMELIST  input  to  the  program.  An  arbitrary  sinusoidal  curve, 

PS0P0I ( time) ,  as  described  near  the  beginning  of  the  program,  is  avail¬ 
able  for  use  as  the  source  of  the  calculation  parameter,  PSOPOI.  This  is 
the  default  form  of  the  PSOPOI (time)  parameter  and  is  a  rough  approxima¬ 
tion  of  the  experimental  periodic  flow.  The  NAME*. ‘ST  input  to  the  program 
must  include  CURVE  =  .FALSE,  if  PSOPOI(time)  data  is  to  be  read  from  a 
data  file.  Input  file  TAPE5  must  include  all  input  data,  both  the 
NAMELIST  variables  and,  if  CURVE  =  .FALSE.,  the  PSOPOI  data. 

The  program  output,  to  output  file  TAPE6,  includes  input  conditions 
and  averages  of  ejector  performance  parameters.  Output  of  performance 
parameter  values  2t  each  time  increment  is  the  default  form  of  the  output 
and  can  be  suppressed  by  including  OUT  =  .FALSE,  in  the  NAMELIST  input. 

The  value  of  MSI  at  which  the  supersonic  regime  breaks-off  with  the 
approximate  regime,  MCRIT,  and  the  corresponding  stagnation  pressure 
ratio,  PCRIT,  are  written  onto  the  output  file  TAPE6  for  each  data  set 
read  into  the  program.  PS1P1  is  written  onto  TAPE7  at  each  time  incre¬ 
ment.  This  feature  was  used  to  obtain  a  theoretical  PS1P1  versus  time 
curve  to  compare  with  experimental  results. 
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Time  is  normalized  by  the  time  increment  between  successive  PS0P0I 
data  points. 
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A. 1  UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 


UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
MAIN  PROGRAM 


PROGRAM  UCAE (INPUT, OUTPUT,  TAPE?  ,TAPE5=  INPUT, TAPB6  =  OOTPOT) 

C  UNSTEADY  CONSTANT  AREA  EJECTOB  INTEGR ATIOH  PROGRAM 
C 

C  THIS  PBOGHAH  CALCULATES  THE  ONE-DIBENSIONAL  T  HE08ECTICA  L  NBA  M  VALUES 
C  0*  THE  EJECTOR/DUCTED  BOLTIPLE  STREABS  PEREOHBA  NCS  PARABETBRS  FOR  THE 
C  GEOBETRtC  AND  PLOB  BOUNDARY  COBDITIOHS  IHPDT.  A  ?SOPO=PSOPO  (TIflE) 

C  CORPS  IS  TO  BE  THE  SOOHCE  0?  THE  CA1COIATIOS  PARABETER,  PSOPO.  THE 
C  CORVE  IS  SUBDIVIDED  IHTO  STEPS  OR  SLOGS  IR  TIRE  AMD  THE  AHALISIS  IS 
C  PEREORBED  OB  EACH  STEP  TO  DETERMINE  THE  AVERAGE  VALUES  OP  THE  PRR- 
C  PORBANCE  PARAMETERS  OP  INTEREST.  THE  PSOPO  (TIME)  CORVE  CAM  BE  IMPOT 
C  VTA  DATA  »TLE  OR  AM  A8BITRART  STROSOIDAL  CORVE  IS  AVAILABLE.  TIBE 
C  IS  NORMALIZED  BT  THE  SAMPLING  TIME  INCREMENT. 

C 

c  ••*•»**•»*••••*••»*•»***«■***»•**«*•*****«•••*•*.••.•,»••• 

C 

C  VARIABLE  DEFINITIONS  AND  THEIR  DEPAOLT  VALDES: 

C  SOME  VARIABLE  DEFINITIONS. 


C 

C  GS,GP....  THE  RATIO  OF  SPECIFIC  HEATS  OF  THE  SECONDARY  AND 
C  PRIMARY  PLOIDS. 

C  BWSP .  THE  SECONDARY  TO  PRIMARY  HOLECOLAR  BEIGHT  RATIO. 

C  TSOPO....  THE  SECONDARY  TO  PRIMARY  STAGNATION  TEMPERATURE  RATIO. 

C  HP  1 .  THE  PRIMARY  PLOW  MACH  AT  STATION  1. 

C  MSI .  THE  SECONDARY  PLOB  MACH  NUMBER  AT  STATION  1. 

C  MM3 .  THE  MIXED  PLOB  BACH  ROBBER  AT  STATION  3. 

C  MBAX .  THE  SAXIBOM  MACH  NOBBER  THAT  THE  PRIMARY  PLOB  CAB  ATT  AIM 

C  BITHIN  THE  NIXING  SECTION  AS  A  RESOLT  OP  A  ONE  DIMENSIONAL 

C  ISENTROPIC  EXPANSION. 

C  PM3S0I...  INPOT  VALOE  0?  THE  BACK  PBESSORB  CONDITION,  RATIO  OP  THE 
C  SITED  PLOB  EXIT  PLANE  STATIC  TO  SECONDARY  STAGNATION 

C  PHESSORE. 

C  PHOSO....  THE  RATIO  0?  SUED  TO  SECONDARY  PLOB  STAGNATION  PRESSORES. 
C  PS1P1....  THE  RATIO  OP  SECONDARY  TO  PRIMARY  PLOB  STATIC  PRESSURES  AT 
C  STATION  1. 

C  PSOPO....  THE  SECONDARY  TO  PRIMARY  STAGNATION  PRESSURE  RATIO. 

C  A  P 1  M3 ...  .  THE  RATIO  OE  THE  PRIMARY  TO  TOTAL  ELOB  AREAS. 

C  MCASE....  A  DATA  SET  COUNTER. 

C 


C  DEPAOLT  VALUES  OP  THE  INPOT  VARIABLES. 

C 

C  GS(1.4)  ,GP(1.4)  ,MBSP(1.0)  ,TS0PO(1.  0)  ,  BP1  ()  ,PH3S0I()  ,  MHA I  ()  ,  AP  1  M3  () 

C 

C  LOGICAL  VARIABLES,  THEIR  PORCTIONS  AND  DEPAOLT  VALUES. 

C 

C  NEB . TRUE.  IS  REQUIRED  EACH  TIME  THE  DATA  SET  HAS  A  DIPFEHEHT 

C  MP1,  API  M3 ,  GS,  OR  GP  THAN  THE  PREVIOUS  SET.  TO  AVOID 


C  CALCULATING  NCRIT  WHEN  HO  CHANGES  OCCUR  IN  THESE  VARIABLES, 

C  5E»=.  FALSE.  SHOULD  BE  INPUT.  .TRUE.  IS  THE  DEFAULT  VALUE. 

C 

C  OOT .  CAOSES  OUTPUT  FOR  EACH  TIBE  STEP  BHEN  .TRUE.  IS  THE  INPUT. 

C  .PALSE.  SHOTTLD  BE  USED  I?  CYCLE  AVERAGES  AND  CASE  INPUT 

C  DATA  ARE  THE  ONLY  OUTPUT  DESIRED. 

C 

C  CURVE....  I?  AN  ARBITRARY  SINUSOIDAL  PSOPO  (T)  CORVE  IS  TO  BE  USED, 

C  .THOE.  MOST  BE  THE  INPUT.  THIS  IS  THE  DEFAULT,  FOR  BHICH: 

C 

C  AMP*. 100,  THE  AMPLITUDE. 

C  P AVG* .100,  THE  MEAN  PSOPO. 

C  HCYC  =  10,  THE  NUMBER  OP  TIBE  STEPS  PER  CICLE. 


0CA001 
UCA002 
0CA003 
UCA004 
UCA005 
UCA006 
UCA007 
UCA008 
0CA009 
OCAOIO 
HCA011 
OCA012 
UCA013 
OCA014 
0CA015 
DC  A  0 1 6 
UCA017 
0CA018 
DC AO  19 
UCA020 
0CAO21 
OCA022 
UCA023 
0CA02U 
UCA025 
0CA026 
0CX027 
0CAO28 
UCA029 
DC  AO  30 
UCA031 
UCA032 
UCA0  33 
OCA034 
UCA035 
UCA036 
0CA037 
UCA038 
DCA039 
UCA040 
UCA041 
UCA042 
0CA043 
UCA9  44 
UCA045 
UCA046 
0CA0  47 
0CAO48 
0CA049 
DCA050 
0CA051 
UCA052 
UCA053 
DCA054 
UCA055 
DCA056 
UCA057 
DCAD58 
UCA059 
OCA060 


on  n  o  on  no  no  n  no  on  non  r.  nnnonnoonnnnnnnnno 
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UNREADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
MAIN  PROGRAM 


ABE  TUB  DEFAULT  VALUES.  OCA 06  I 
If  A  PSOPO  (T)  CUBVE  IS  TO  BE  BBAO  flOH  A  DATA  fits.  . f  ALS A.  UCA062 
BUST  SB  IBPUX.  UCAOoJ 

OCA  064 

OUTPUT  DATA:  UCA065 

OCA 066 

AU  DATA  IMPUT.  CICLE  SAAB  VALUES  Of  BS I ,  PSlPl,  HH3,  PSOPO.  PHOSO,  JCA067 
ABD  BSP.  BOTE  THAT  OUT-. TBUE.  BILL  CAUSE  TUB  ABOVE  TO  BE  OUTPUT  PUB  UCA06B 
EACH  TIUB  STEP  II  TBS  PSOPO(TIHE)  CICLB  AS  BELL  AS  THE  BBBAK-OPP  UCA069 
VALUES  OP  PH3SO,  PHOSO,  AID  HH 3  IB  SB,  APPU,  3 SB  PLOWS,  ABO  HB  UCA070 
PLOW  SOLUTIOBS  THAT  BEET  TUB  SACU  PBESSUHE  COBOIIIOB.  UCA071 

UCA 072 
UCA073 
UCA 0  74 
UCA  075 

HUC U  OP  THIS  PBOGBAH  COSES  OIBECILI  FBOB  CAE,  COISIABT  ABEA  EJECT OB  UCAQ76 
PBOGBAH.  USED  IB  "AU  IUVESTIGAXIOU  OP  SAS-OIUA&IC  PLOW  PBOBLEHS  IB  UCA077 
CUEHICAL  LASEB  SISIEHS",  AUDI,  A.  L.  ABO  BIKKELSOB,  C  D. ,  UCA078 
UBI VEBSII X  OP  ILLH01S  BEPOBI  UILU-EBG-74-4009,  1974.  UCAJ79 


UCA 080 
UCA08 1 
UCA082 
UCA08J 
UCA064 


COHHOB/BLOCE 1/GP, GS, HP  1 , HUSP,ZSOPO, AP IHJ.EHIB , UEB, HHAI ,  UCA OSS 

SPCBIT  UCA086 

COHHOB/BLOCK2/HSl,afl3, PS 1?1. PHOSO, BSP, PB3SO  UCA087 

COHHO B /BLOC  X3/5H5 1 . SPHOSO. SHHJ.SPSOPO . SPS I P 1 , SB  SP  UC  AO  UH 

LOGICAL  CUB VE, BEU, OUT  UCAUd9 

BEAL  US  1  ,H P  1,HCBIT, SUSP,  H,  SB 3,  BUB,  H BAX  UCA090 

0CA091 
UCA  092 

HABELIST/I  UCAE/GS,  GP,BBSP, TSOPO.AP  1H3,  HP  1,  PH3S0I  .OUT , BCASE  UCA093 

U.HEU.SHAX, CUBV2, AHP, PA VG, BCIC  UCA094 

ASSIGB  DEFAULT  VALUES.  UCA095 

DATA  GS.GP, SUSP, TSOPO, OUT, BE B, BCASE/ 2*1. 4,2* 1. , 2*. TBUE. , 1/  UCA096 

DATA  CUBVE, AHP, PAVG, BCIC/. TBUE. ,0. 1U , 0. 10,10/  UCA097 

BOTE  THAI  BCASE  IS  BOX  IBCBEBEBTSD  81  THE  PBOGBAH.  UCA098 

UCA  099 

PPO(G.H)-  <1.0*.5*<G-1.Q)  "(H*fl)  )*•  <-G/  (G-1.  ))  OCA  100 

105  IPHESS-0  UCA 101 

THE  VARIABLE,  IPBSSS,  IS  THE  TIBS  STEP  COOBTEB  IB  THE  PBOGBAH.  OCA  102 

I SPOT  PLOW  ABO  GEOHETIIC  COBDITIOBS.  UCA103 

BE  AO  (5  ,IUC  AE)  UCA104 

IP  (EOF  (5) )  999,110  UCA  105 

110  If(.BOT.BSB)  GO  TO  115  UCA  106 

PBOGBAH  SBBCS  IBS  BAZIHUE  BS1  WITH  BUICH  TUB  PABBI  CHOKING  COBDITIOB  UCA  107 

CAB  BE  SATISFIED.  IBIS  IS  A  PUBC3I0B  OP  BP1,AP1B3,  ABD  GP.  UCA108 

CALL  CBITH  (HCBIT)  OCA  109 

TBS  PSOPO  COBB  ESP  OBOISG  TO  HCBIT  IS:  UCA  110 

PCBIT*PPO(GP,HP1)/PPO(GS,aCBIl)  UC A  11 1 

TUB  IHPUT  IS  BOH  OUTPUT  UCA112 

115  BBITB(6,  120)  BCASE, GS,GP,HUSP,TSOPO, API a3,BP1, PH 3SOI  UCA113 

120  POBHAT (5 (/) ,10X,"IBP01  DATA  FOB  CASE  BO. ", 13 ,//, 101, "GS=", F7. 4 ,  UCA114 

6  8X,"GP-",P7.4,dI,"HUSP-",P7.  4,//.  10X.  “TSOPO-"  ,  F7.4 , 51,  "AP  IB  J«"  ,  UCA  115 

6P7.4.5X,"HP  1»",P7.4,//.  10  A"PH  JSOI"" ,  PS.  4  ,//)  UCA  116 

IP  (N SB)  UBITE  (6,125)  HCBIT, PCBIT  UCA  1 17 

125  POBHAT  I//.  10X,  “FOB  THIS  CASE,  HCBIT- " ,P7.4 ,5X. "PC8IT-" , f 7.4, //)  UCA118 

UCA)  19 
UCA120 


A 
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UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
MAIN  PROGRAM 


130  CALL  PR AT TO (ASP, PA  VG , NCTC, PSCPQI, I  PRES  S , COR  VE)  0CA121 

IF(IPRSS?.FQ.NCYC.  AND.  •  ROT.  COP  V  g)  GOTO  145  OCA122 

IF (.NOT.CORVE.OR.IPRESS.GT.1)  GO  TO  140  7CA123 

WRITS (6 , 135)  ABP.PAVG, HCYC  0CA124 

135  FORBAT (/, 10T,"FOR  THE  SIROSOIDAL  PR  ESSO  RE  RATIO  C3RVE; 10 X,  OCA125 

G"ASP*",  "7.  3,  10  X,  "PAFG*",F7.  4 , 1 0  X,  "NCYC*",I5,  "  TIB?  SLQGS/C  YC  LE")  9CA126 
140  CALL  CALCAE (PSOPOI.PNISOI.OOT.IPBESS.NCYC, NEBROR.CORYE, SCR IT)  0CA127 

I»(NSRROR.GT.O)  GO  TO  105  0CA128 

IF  (IP8ESS.LT.NCYC)  GO  TO  130  0CA129 

145  WRIT5(S,  150)  NCTC  OCA130 

150  FORSAT (//.10T.-THE  TOTAL  ROBBER  OF  TIHE  SLOGS  BEAD  WAS", 15  ,////)  OCA131 

GO  TO  105  OCA132 

494  CONTI  ROE  0CA133 

END  OCA134 


■ion 
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UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  CALCAE 


CALOOI 
CAL002 
CAL003 

c  CALOOtt 

SUBROU'INE  CALCAE(PSOPOI,?H3SOl.OOT,IPRESS,NCYC, TERROR, CORTE,  CAL005 

SHCEIT)  CAL006 

C  THIS  SUBROUTINE  DETERMINES  THE  “LOW  REGIME,  CALLS  THE  PROPER  SOB-  CAL007 

;  “0  "TINES  TO  DETERMINE  THE  PLOW  CONDITIONS  AT  STATIOKS  1  C  3,  AND  CAL008 

C  CONTROLS  ITERATION  PROCEDURES.  CAL009 

-  CAL010 

C  CAL011 

COMMON/BLOCK  1 /G P. GS, N P 1, H »SP. TSOPO, API 13," MIN, M MAX, NEW,  CAL012 

6PC“IT  CAL013 

COMNON/3LOCK2/HS1,HM3, PS1P1, PHOSC, WSP.PM3SO  CAL014 

“ESL  '*S1.*P1,MCRIT,NWS?,N,Kr,NM3,MS1B,MR,MMIN,MNAT  CAL015 

LOGICAL  NEW , CU  R V  E, OUT  CAL016 

-  CAL017 

C  w*  SQRT  (“*?0)  /?“A“WN  (G,H)  CAL018 

WN  {1,  HI  =N*SQHT  (G-  (  1.  3*  (G-  1.  7)  *  (N*M)  *.50)  )  CALO 19 

C  TH“  P/PO  ?"MCTION.  CAL020 

PPO  (G  ,  N)  =  (1 . 0*  (G-1. 0)  *  (N«S)  *  .  50)**  (-G/  (G-1. 0)  )  CAL021 

7  THE  INVERSE  OF  THE  P/PO  “UNCTION  CAL022 

PH  (G,  PS)  “SQRT  ( (PR**  ( ( 1. 0-0)  /G)  -  1. 3)  *2.  0/  (G-1 . 0)  )  CAL023 

^  ot/oy  *or  !  NORMAL  SHOCK  W V 7  E .  CAL024 

PYYMX  (G,.NY)  =  (2.0«G/(3»1.0)  )  *  (MX*HT)  -(  (G- 1 . 0)  /  (G  *  1 . 0)  )  CAL025 

CAL026 

c  CAL027 

DA  ?A  ERR.NS1,  ERHORT,  NITHAX/2*  1.  OE-4,  35/  CAL028 

OAT'-  SR/”SP,,/,SSa/,,SSR"/»NR/"H5  V. APPP/"APPR"/  CAL029 

C  MS  IB  IS  THE  LOWEST  SECONDARY  EACH  ROBBER  CONSIDERED  IN  A  SEARCH  FOR  CAL030 

-  1  SOLUTION.  CAL031 

N“RROR=0  t  NCOUNT=1  S  DMS1B“.100  $  NS1B*.001  CAL032 

C  ~H“  NAYI.NUH  PS1P1  ALLOWED  “OR  CALCOLATIONS  TO  PROCFDE,  PS1P1U,  CAL033 

C  IS  BASED  ON  A  NORHAL  SHOCK  STANDING  IS  THE  EXIT  PLANE  OF  THE  CAL034 

C  PRI.NART  NOZZLE.  CAL035 

PS1°1r'““Y*'NT  (GP,“P1)  CAL036 

AS  1  P  1  “  (1.0/AP1Hl)-1.0  CAL037 

-TNST  =  SORT  (NWSP/TSOPO)  *  AS  1  P1/WB  <GP,MP1)  CAL038 

C  I"  THE  RS3INE  IS  SSR  PS1P1>“  1.0  AND  BS1=  1.0  CAL039 

c  CAL040 

RFGI’,2“SSR  CAL041 

US  1  —  1 • t  CALO  4  2 

PS1P1  =  ?PO(GS,  1. 0)  *PSOPOI/PPO  (GP,HP1)  CAL043 

I" (PS1P1.LT. 1.0. AND. PSOPOI.LT. PCRIT)  GO  TO  20  CAL0U4 

IF(?S1?1.GT. 1.0)  GO  TO  50  CAL045 

C  THE  PRESSURE  RATIO  LIES  BETWEEN  SR  AND  SSR  REGIMES.  PS1P1  =  1.0  CAL046 

C  IS  A  GOOD  APPROXIMATION  CAL047 

REGT“E“  A  ?°R  CAL048 

JS1?1“1.0  CAL049 

PS1?O=1.0*PPO(GP,MP1)  /PSOPOI  CAL050 

MS1=PN  (GS,  PS  ISO)  CAL051 

GO  -0  50  CAL052 

C  IN  SR  ESI  MUST  BE  WINDOWED  SUCH  THAT  NS  1 8  <  *  NS  1  <  NCRIT  CAL053 

20  REGINE  =  SR  5  NI~=1  S  NTT  PE=  1  CAL054 

MS1»NS1B  CALO 55 

10  CALL  “ABRI  (MSI  ,  PS1P1 ,  NERROR)  CAL056 

T“  (NER’OE.  GT.  0)  GO  TO  906  CAL057 

PS1P1 I  =  PSO?OI*PPO(GS,  MSH  /PPO  (GP,  SP1)  CALO 58 

I“(  ( NI“ . “Q.  1)  .AND.  (PS1P1.GT.  PS  1  PIT) 1  GO  TO  909  CAL059 

Y*  (PS1P1-PS1P1I1  /PS1  PIT  CAL060 
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UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  CALCAE 


IP  1  (ABS  (X)  .GT.  BBBOBI)  .  AMD.  (PS  IP  1.  LB.  PS1P1I) . AMD. 
s  (ass (asi-dtaiT)  /acaii.Li..  ooo5) )  go  to  40 
IP  (  (SSI ♦DBS  18)  .GT.BCBIT)  DBS1B»  (BCBIT-HS  1)  /2.  0 
CALL  IX£B(aS1,DaSlH.EBBBSt.+  1.  0,1.0.  o.ebbobx.bi i.bttpe, 

6XNSG,  YNSG.XPOS.XPOS  ,  NSIGM 1 ,  M  SIGB2) 
IF{HirflAX.LI.NIT.ABD.MTIPE.KE.3)  GO  TO  40 
GOTO  (30,30,50)  .STIPE 

C  IP  BO  SOLUTION  HAS  SSLS  POUND  SB  BILL  BE  SEABCBEO. 

40  IP(OUI)  MBITS (6, 906) HS1,PS1Pl,PS1P1I 
GO  TO  SO 

50  CALL  OCV(HEBBOB) 

IP  (BBBBOB.BE.  0)  GO  TO  902 

C  IP  THE  VALUE  OP  PB3SO  BETUBNED  SI  OCV  IS  <  PH3SOI  THE 
C  SECONDAHI  BASS  PLOBBATE  IS  DEPEBOEBT  OB  PB3SOI  ABD  TEE  BEGIBE 
C  BUST  BE  H  B. 

DELT=*PB3SO-PB3SOI 

IP  (BEGIBE.  BE.  BB.  ABO.  PB3S0. GT.PB3SOI)  GO  TO  70 
IP  (ASS  (DEL  T/PB3SOI)  .  LI.  EBB  BS  t.  ABD.  BEGIBE.  aQ.  BB.  AND. 

SDELT.  LE.  (0.0))  GO  TO  70 

THE  BACK  PBESSUBE  IS  TOO  BIGH  FOB  THE  SECONDABX  TO  BE  CHOKED. 

BS 1  IS  LESS  THAN  THE  BBEAK-OPP  VALUE  ABD  THE  BEGIBE  IS  BB.  A 
VALUE  OP  BS  1  PHODUCING  A  Pfl3SO» PB3SOa  IB  OCV  IS  SOUGHT. 

60  IP  (BEG IBs. B E. BB)  BS  1^.500 
BEGIBE-BB 

IP  ( BCOUBT.  GT.  25.0B.  BS  1.LI.  BS  1 B. OB. BS  1.GI..  999)  GO  TO  904 
CALL  I IB Si (HS1.PH3SO.PH3SOX, BCOUBT, BS 1 LO, SS 1HI) 

PS1P1»PSOPOI*PPO  (GS,BS1)/PPO  (GP.BP1) 

GO  TO  50 

70  IP(PSIPIU.LT.PSIPI)  GO  TO  900 
HSP»CONSI*PS1P  1*BH  (GS,  BS1) 

IP  (.NOT. OUT)  GO  TO  100 
IP  (IPEESS.  EQ«  1 )  HBIIE  (6,00) 

80  FOB SAT  (//, 10X, "TIBS  SLUG  DAT A",/, 3X, "SLUG", 2X, "BEGIBE", 4X, 

5  "PSOPO  ",6X , "BSP" ,6  X, "BS 1", 5X,"PS1P1*»5X,"flB3",5I, 
6,"PBOSO".7X."PH3SO"./) 

BBITE  (6,90)  IPHESS.  BEGIBE, PSOPOI, BSP  ,  BS1 , PS  IP  1.HH3 

6  ,  PflOSQ,  PB3  SO 

90  POBBAT  (/, 2X.I4, 5X,  A4.3X ,P6. 4.4X.P6.4, 3(31,  P6.  4)  , 

&3X.P7.  4.3X.P8. 4) 

BBITE (7.95)  PS  IP  1 
95  POBBAT  (10P10.4) 

100  CALL  SUBUP  (IPEESS, HCXC, PSOPOI, CUBVB) 

BE  TUBS 

c  ••*•**••«  ESBOB  BESSAGES  **••**»*•* 

900  MBITS  (6,901)  IPBESS,  PS  IP  1.  PS  IP  10,  BS  1 

901  POBBAT  (/."•  ***•••*•«•***•*•  ",/,"EBBOB  IB  CALCAE 
S"./,"*  •***•*  ***••**•  ",/,5X,"TIBB  SLUG  NO.",I4,5X," 
&PS1P  1*"»P8.  4,5X,*PS1P1U*"»P8.4,5X,*BS1"",F8.  4) 

BEBBOB»1  t  BETUBB 

902  BBITE  (6,903)  BS1.PS1P1 

903  POHBAT (//,"*  **»**••***••  »,/,«KBBOB  IB  OCV",/,"  * 

6*  *••*•*•  *«./,10X,"HSl*",F8.4, 1QX,"PS1P1»",P8.4) 

H8TUBM 

904  HBITB  (6,905)  BS  I . PB3SO,  PB3 SOI 

905  POBBAT (//."*  ************  *",/,«BBB0B  IB  CALCAE",/,"  • 
&•••••••  ",/, "MOM—CQBV EBGEBCE  OP  BS 1  PUB  THE  INPUT  PH JSO," 

S,1X,"HB",/.10X,"BS  1*",P8.4,  10X,  "Pfl3S0»", P8.  4 , 10X, "PB3S0I"",  F8.  4) 

BEBBOB* 1  S  BETUBB 


CAL061 
CAL062 
CALD63 
CAL064 
CAL065 
CAL066 
CAL067 
CAL368 
CALD69 
CAL070 
CAL071 
CALU72 
CAL073 
CAL074 
CAL075 
CALO  76 
CAL077 
CAL078 
CALO 79 
CALJ80 
CAL081 
CAL 082 
CAL083 
CALO  84 
CALQ85 
CALO 86 
CAL087 
CAL 088 
CAL089 
CAL090 
CAL 091 
CALO 92 
CAL093 
CAL094 
CAL095 
CAL096 
CAL 097 
CAL098 
CAL099 
CAL  100 
CAL  101 
CAL  102 
CAL  103 
CAL  104 
CAL  1 05 
CAL  106 
CAL  107 
CAL  108 
CAL  109 
CAL1 10 
CAL1 11 
CAL  112 
CAL1 13 
CAL1 14 
CAL  1 15 
CAL  116 
CAL117 
CAL  1 18 
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UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  CALCAE 


°0S  JPI ■'?  ( 6 , 907)  CAL11B 

q07  -r)?-M.T  (/,"  -  *»,/, "ERROR  r:i  7ABRI",/,  "  •  *  •  CAL  1 1 9 

5  *•*•••  •  *")  CALI 20 

M?RP9R=’  $  RFTSRN  CAL121 

q  iq  •»'>•>"*  :  (///,  13X,  "SOLTTICM  NOT  PC- OS!)  IS  SR  DY  MIT!!  IT2P",/,  CM122 

5l''<,',SSl  =  ",?3.5,h:<("?Slpl^»,  ?S.  5,  6X,  "PS  1?  1 1*  •  3. 5)  CALI  23 

009  BRIT'’ (6,913)  CAL124 

010  ’o’-r  p  (/I  ,  5X,  "PSOPOI  IS  SO  LOM  TH  S.T  "SK.331  IS  SR.")  C1L125 

!lr?->OR=l  5  RETOSS  CAL126 

I”9  CAL127 
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UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  FABRI 


c  **•.•*•••****•••**•*••«•*•••***••**•*•••*•*•••***•••*•••■*.* 

r* 

c 

subroutine  fabri <hsi,  psipi.nehbori 

c 

C  CONSTANT-AREA  EJECTOR ,  FABRI  CRITERION  (FABBI)  SUBROUTINE. 

n 

C  1-D  EJECTOB  ANALYSIS  WITH  SUPERSONIC  PRIHAHT  (HP1>1)  AND 
C  SUBSONIC  SECONDARY  (HS1C1)  WHICH  CHOKES  IN  THE  EJECTOR 
C  SHROUD  (*S2=1) •  THIS  EJECTOB  ANALYSIS  IS  BASED 
C  ON  THE  CONS~ANT-ABEA  SHROUD  EJECTOB  ANALYSIS  OF  FABRI. 

C  REFERENCE:  HACA  TN  1410. 

C 

TOSH  0*1 /BLOCK  1/GP  ,SS,NP1#HWS?,TSOPO,APlH3,NRIM,NEW.HR\Xf 
SPC3IT 

REAL  H.HP1, NS1, HWSP, HH IN , NH A  X , N P2 
LOGICAL  NEW 
C 

C  ?/PO(J,Nl  FUNCTION. 

PP0(G,H)*(1.f.5«  (G-1.)  )  **  (-G/IG-1.)) 

C  FUNCTION. 

AA3H(G,N)*  (1./B)  *(  ((2./(G*1.)  )  *  (1 . *0 . 5*  (G-1  .)  •  (H*N) )  ) 

S  **(0.  5*(G»1.)/(G-1.))) 

C 

C  FABRICS  CRITEBION  APPLIES  ONLT  WHEN  (BS 1.  LB.  ( 1.  0)  )  AND 
C  (PSIP’.LE.  (1.01  )  . 

C  CALCULATE  THE  AREA  RATIO  FOR  STREAH  (P)  AT  STATION  (2). 

AP  2PS  *  (  A  ASH  (GP.NP1)  /AP1N3)  •  (  1  .  -  (1 .  -  API  B3) /A  ASH  (GS.HS1)  ) 

C  CALCULATE  NP2  AT  STATION  (2)  WHERE  HS2  =  1.0. 

CALL  ASTAR  (GP.  AP2PS,  HP2  ,  HHIN,  HHAX) 

C  CALCULATE  ?“A~IC  PRESSURE  RATIO  PS  IP  1  BASED  ON 
C  FABRI"  S  C?  IT  SR  ION. 

YD=  (1. 0-AP1H3)  •  (  ( 1. 0*GS«  (NS  1**2)}  - 
S  (PPO  (GS.  1.0)  /PPO  (GS.NS1)  1  *  (  (1.0  +  GS1/AASH  (GS.HS1)  )  ) 

VN*  (PPO  (GP,  HP2)  /PPO  (GP,  HP1)  )  «AP  IN  1*  (A ASH  (GP,HP2)  / 
r,  AASN  (GP,  HP1)  )  «  (1.  *GP*  (NP2«*2))  -  AP1H3*  (1.  *GP*  (HP1*HP1)) 
PS  1 P 1*7  N/VO 

I"(PS1P1.GT.  1.0.  AND.  .  HOT.  NEW)  WRITE  (6,80) 

IF(PS1P1.LE.  (0.0)1  WRITE  (6,6  0) 

I  ®  (HP2  .  LT.  N  P 1)  W  RITE  (6 , 40) 

I?(PS1?1 .LE. 0.0. OH. HP2.LT. API)  GO  TO  90 
I'tPSIPl.G".  1. 0.  AND.  .NOT.  NEW)  GO  *0  90 
RETURN 


c 

ERROR  HESS  AGES. 

40 

•ORHAT ( / , 5 X , " .  . 

.  ERROR 

IN 

F  ABHI : 

(SP2.  LT.  HP  1 )  .  ",//) 

60 

•ORHAT (/,5X , " • . 

. ERROR 

IN 

FABRI: 

(PS1P1.LT.  (0.0)  )  ." 

8" 

FORHAT (/,5X, 

.  ERROR 

IN 

•ABHI: 

PS1P1>  1.0",//) 

c 

SET 

F8ROR  »LAG. 

90  NEPROR=2 
RETURN 
END 


FA8001 
FAB002 
FAB003 
FAB004 
FAB005 
FAB006 
FAB007 
FAB008 
FA BO 09 
•ABO  10 
FAB011 
FAB012 
FAB013 
FAB014 
FAB015 
•AB0 16 
FAB017 
FAB018 
F AB0 19 
FAB020 
FAB021 
•ABO  22 
FAB023 
?  ABO  2  4 
•AB025 
FAB026 
FAB027 
FAB028 
FAB029 
PAB030 
FAB031 
•ABO  32 
FAB033 
FAB034 
FAB035 
•AB036 
FAB037 
•AB038 
•ABO  39 
FAB040 
FAB041 
P‘B042 
•AB043 

FAB044 
FAB045 
FAB046 
PAB0  47 
FAB048 
FAB049 
FAB050 
FAB051 


Sm-'P1*c£.  > 
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UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  OCV 


C  ••••••••••••••••••••••••••••••••••••••••••»••**•••••••••••**’ 

c 

SUBROUTINE  OCV(NERBOR) 

C  CONSTANT  AREA  EJECTOR  OVERALL  CONTROL  VOLUME  ANALYSIS  (OCV)  . 

C  1-0  ANALYSIS  'ROM  INLET  STATION,  1,  TO  ETIT  STATION,  3. 

C 

C0MM0N/BL0<'K1/GP,GS,HP1.MVSP,TS0P0,AP1H3,HMIM,NEW,HMAX, 
SPCR IT 

COM MON/BLOCK  2/MS  1 , MM  3 , PS  1 P 1 , PRO  SO , NSP , PH 3S0 
REAL  N  ,  M  "  1 ,MS1 , MM  3 , MM 3M , MM 3P , MMIH , M WSP , HWM P, HSQD 3H . HSQ D3P 
6,HMAX,HR 
LOGICAL  NEW 

WM  (G  ,  M)  «M*SQHT(G*  (1.0*0. 5*  (G-1.0)  *  (H*N)  )  ) 

PPO(G.M)  =  (  1. 3*.  5*(G-1. 0)  *  (M*M)  )  **  (-G/ (G-1.0)  ) 

C  A  FORM  0”  THE  1-D  THRUST  'UNCTION. 

T  (G,M)  =  ( 1.  3*G*  (3* MM  /(M'SOBT  ( 1 .  0+ .  5»  (G-  1 . 0)  •  (M* H) ) ) 

C  REDEEINE  MR  LOCALLY. 

.MR  =  ',3R" 

C  CALCULATE  CONSTANTS. 

CD=SORT  (MWSP/TSOPO) 

C  'J  ECTOR  MASS-'LOW  RATIO 

BSP*P«  1  P1»  (  (  1.  -  AP1M3)  /  AP1M3)  *C0'{WH  (GS  ,  MS  1 )  /W H  (G P,  MP1)  ) 
CPSP  =  (G3/GP)  *  (  (GP-1.)  /(GS-1.)  )  /HHSP 
C  MIXED  'LOW  PROPERTIES 

MWMP»  (1  .  *WSP)/(1  .♦  (WSP/MWSP)  ) 

0M»1.  /  (1.-  (  (GP-  1.  )  /GP)  *  (  (1.  ♦  (WSP/MWSP)  )/  (1.  +CPSP*WSP)  )  ) 

Z  CALCULATE  THE  MITED-TO-P5IMARY  STAGNATION  TEMPERATURE  RATIO. 

TMOPO=  (  1 .  ♦  WSP'CPSP'TSOPO)  /  ( 1  . *WSP*CPSP) 

Z  CALCULATE  SOME  CONSTANTS. 

C1*SORT  (  (TSO PO/MWSP)  *  (GP/GS)  ) 

~2  »SQR'"  ( '"30  PO/3  WM  P*  (GP/GM)  ) 

C  SOLVING  'OR  M33. 

-Ml*  (T(GP,MP1)  ♦C1*WSP*T(GS,MS1)  )/{  (l.'WSP)  «C2) 

TM1MIN  =  SDRT  (2.*  (GM'1  .  )  ) 

IF  (TM3.  LT.  T33MIN)  GO  TO  20 
C3’(T',3**2-2.*GM) 

C4=  (  (  (GM-1.)  «.5)  •  (TM3*«2)  -GM*GH) 

C5’e,0R”(C3*C3+U.  *C4) 

3S3D3M=  (-C3-C5)  /  (2.«C4) 

MSJD3P=  (-C)*C5)  /  (2.  *C4) 

C  DETERMINE  TWO  POSSIBLE  MIXED-FLOW  MACH  NO.  SOLUTIONS. 
I'(3S0D3H.GE.  (0.0))  MM  3M«SQ8  T  (8  SQD3H) 

I'(MS0D3P.GE.  (0.0))  MM3  P'SQRT  (H  SQ  D3P) 

C  USE  ONLY  3U3S0NIC  RESULT  AT  (3)  . 

MM3=MM3P 

C  CALCULATE  STATIC  PRESSURE  RATIO,  PM3P1,  AT  (3). 

CM  =  SORT  (  TMO  PO/M  WMP) 

PM3P1»C6*.AP1M3*  (1.  *WSP)  •  (WM(GP,  MP1)/WM  (GM,  MM3)  ) 

C  CALCULATE  OTHER  STATIC  AND  STAGNATION  PRESSURE  RATIOS . 

?M3SO=  (pun  (GS,MS1)  )  *PM3?  1/PS  1P1 
PMOSO=PM3SQ/PPO (GH, MM3 ) 

RETURN 

C  ERROR  MESSAGES. 

AO  WRI'E  (6,30) 

3"  'ORM’T  (/,5ufn.  ..  ERROR  IN  OCV:  NC  SOLUTION.",/) 

NET  ROR= 1 

RETURN 

END 


ocvooi 

OCV002 

0CVD03 

OCV004 

OCV005 

OCV006 

OCV007 

OCV008 

OCV009 

OCVOIO 

OCV011 

OCV012 

OCV013 

OCV014 

OCV015 

OCVO 16 

OCVO  17 

OCV018 

OCVO 19 

OCV020 

OCV021 

OCV022 

OCV023 

OCV024 

OCVO 25 

OCV026 

OCV027 

OCV02B 

OCV029 

OCVO 30 

OCV031 

OCVO  32 

OCV033 

OCVO 34 

OCV035 

OCVO 36 

OCVO 37 

OCVO 3 8 

OCVO 39 

OCV040 

OCV041 

OCVO 4 2 

OCV043 

OCV044 

OCV045 

OCV  046 

OCVO* 7 

OCVO 4 8 

OCV049 

OCVO  50 

OCV051 

OCV052 

OCVO 53 

OCVO  54 

OCVO 5 5 

OCVO 56 

OCVO 57 

OCVO  58 

OCV059 

OCV060 


non 


133 


UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  PRATIO 


•**•«*•••«••*••••••••••••*•*••*•••••••*••••••••••••••••••••••*  PBAOOI 

PBA002 
PBA003 

SUBBOUTIBB  PBATIO  ( A  HP.PA  VG  ,  BCIC,  PSOPOI.IPBBSS.COBV  E)  PBA004 

C  TUX S  SOBBOOTXBE  IBITIAIIZES  TUB  IBPOT  STAG BAIIOB  PfieSSOBE  BAXXO,  PKA005 

C  PSOPOI.  FOB  BACH  IX HB  SLUG.  PBA006 

C  PBA007 

coauou/BLOct2/as i.auj.ps ip i.paoso.gsp, ph Jso  pbauob 

COHBO  K /BLOCK3/SHS 1 , SPHOSO. SHHJ.SPSOPO.SPSIP 1.SBSP  PBA009 

LOGICAL  COB  VS  P8A010 

PI-3.  141S9J  PBA011 

XP(COBVB)  GO  TO  80  PBA012 

BEAD (3.50)  PSOPOI  PBA013 

IP«OF  (5))  70,60  PB A0 14 

50  POSH  A  I  (F10.  4)  PB  AO  15 

60  IPBBSS— IPBESS+ 1  PBA016 

BC IC-IPBBSS  +  1  PB  AO  17 

BBTOBB  PBA018 

70  BCYC-IPHSSS  PB  A0 19 

CALL  SUHUP  (IPBBSS.BCIC. PSOPOI, C0B7K)  PHA020 

BBIOBB  PBA021 

SO  IPB£SS-IPBBSS»1  PHA022 

PBESS-  FLOAT  (IPBESS)  PBA023 

CXC-FLOAT (BCTC)  PBA024 

PSOPOI »PAVU*{1.0*AHP*SIB  (2.  0*PI*PBBSS/CIC)  )  PBA025 

BBTOBB  PB  AO  26 

BBO  PBA027 


non 
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UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  CRITM 


SOBROUTISE  CRIT8  (HCRIT)  _ 

'BIS  ST’POO'INE  IS  USED  TO  ®IWD  THE  SECONDARY  BACH  NOH9ER,  BS1, 
AT  WHICH  PS1P1  *  1.0  IS  SRCOUNTERED  I*  THE  SH  SOLUTION. 

CO®BOH/3T,OCK1 /GP  , GS , B  P 1  .BWSP.TSCPO#  API  B3  , flfllH  ,N EB, HBAI , 

RPC SIT 

®EU  B?1.aCRIT.BHSP,BHIM,BHAT.B 
LOGICAL  HEW 
** .001 
D  B*  .  1 1  0 
S^YPE* 1 
DO  10  1=1,300 
SBIH=BP1 

TILL  ®AB9I fl, PS1 PI .HERROR) 
lE(PSIPI.GT.I.O)  GO  TO  9 
I®  (H.GS.  ( 1. 0)  )  GO  TO  11 
OLDH*B 
GO  "0  10 

9  I®  (STYPE.  7Q  .2)  00  TO  11 
OH*.  GOODS 
IsOLD' 

NT YPE=  2 
10  B  =  B  ♦  D  H 

n  bc»it=  (a*oLDa)/2.o 

i® (bcrI'.ge.  (i. 0) )  aaiH=api 

P®Tn®E 

END 


CRH001 

CRH002 

CRH003 

CRH004 

CRB005 

CRB006 

CRB007 

CRB008 

CRB009 

CRB010 

CRB011 

CRH012 

CRBO  1 3 

CRHO  14 

CRBO  15 

CRH016 

CRB0 1 7 

CRHO 18 

CRB019 

CRH020 

CRHO 21 

CRB022 

CRH023 

CRH02# 

CRHO 25 

CRH026 

CBH027 

C3H028 

CRB029 

CRB030 


UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  ASTAR 


-  ASTOOI 

c  AST002 

-  AST003 

SUBROUTINE  A  S’- A  R  (G .  A  BA  TIO,  HP2 , 3  S3  H.  SHA  X)  AST004 

C  “HIS  SaBROtJ“TNS  SOLVES  TOR  HP2  GIVER  AP2/AP*  EROH  PABBI.  AST005 

S3.\L  HP2,HHI,HHI!I.H.I1LO,  *1HAT  AST006 

C  “HE  A/A*  ’•tJSC“IO»  A ST007 

A  ASH  (G  #  HI  *  ( 1  •  O/H)  "  (2  •  0*  (1 . 0* .  5*  (G-1 . 0)  *  (H-H) ) /(G+1.  3)  )  ASTG08 

5**  ( (G*  1. 0)  /  (2.0*  (G-1. 91 )  1  AST009 

0  "AKE  SURE  THAT  RUST  IS  BIG  5BOUGH  ?OR  THE  DATA  SET  IHPDT.  AST010 

HHI*HHAT  A37011 

Ht9»HHT'l  AST012 

HIT*1  AST013 

10  AHt  =  AASH(G,  HHI)  AST01# 

ALO-AASH  (G.HLO)  AST015 

f1P2=HL0  +  (ARATIO-ALO)  «  (HHI-HLO)  /  (AHI-ALO)  AST016 

I”  (HIT.  GT.  30)  RETHRH  AST017 

DI?r=*ARATIO-AASH  (G  ,  HP 21  AST018 

IE(>  »«  (DIE*)  /ARATIO.  LB.  .  00001)  RETURR  AST019 

X?  (DIFE.GT.  0.0)  GO  TO  20  AST020 

HHt=HP2  AST021 

GO  TO  30  AST022 

00  HLQ- HP2  AST023 

30  NIT=HI7»1  AST024 

GO  TO  10  AST025 

2ND  AST026 


non 


UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  ITER 


•  XT  BOO  1 

ITB002 
XTB003 

SUBBOUTINE  ITEB  (X,  DX  .EBBOBX  ,SI  GN  ,  X  ,  XGIVEN  ,  EBBOBX  ,  NI  X  ,  NTYPE,  IIB004 

6XMEG.INEG.IP03,IP0S,MSIGN1, NSIGN2)  ITBOOS 

C  SUBBOUTINE  PeBFOBHS  AM  IT2BATI0N  TO  FIND  X  SUCU  THAT  TUB  J.IB006 

C  ABSOLUT*  VALUE  OF  (I- XGIVEN)  IS  .LE.  TO  EBBOBX  OB  TUL  ITB0Q7 

C  ABSOLUTE  VALUE  OF  (X  (I*1)-X  (I)  )  IS  .LE.  TO  EBBOBX.  ITB008 

C  VAMIABLES:  XTB009 

C  X  *  INDEPENDENT  VABIABLE  ITBOIO 

C  OX  -  INCBiatNT  IB  INDEPENDENT  VABIABLE  XT  B0 1 1 

C  EBBOBX  »  aAXIBOa  VALUE  OF  ABS  (X  {!♦  1 ) -X  (I) )  FOB  SOLE  ITBJ12 

C  SIGN  »  *1.0  OB  -1.0,  */~  INCBEH EBT1BG  FBOH  XBUIAL  X  ITBJ13 

Cl  *  DEPENDENT  VABIABLE  XTB014 

C  EBBOBX  »  XAXiaUU  VALUE  OF  ABS  (X-XGIVEN)  ITB01S 

C  XGIVEN  *  GIVEN  VALUE  OF  DEPENDENT  VABIABLE  IIB016 

C  HIT  *  INCBE8ENT  NUBBEB  ITB017 

C  NTXPS  »  1 — INCBiBENT,  2 — XNTEBPOLATION,  3 — SOLUTION  ITB018 

C  IT  B0 19 

DX3*!— XGIVEN  ITBOIO 

IF  ( (ABS  (  DX)  -EBBOBX)  .LE.  0.0)  GO  TO  90  ITB021 

XF(DX.GT.O.O)  GO  TO  20  ITB022 

NSIGN2*— 1  ITB023 

XNBG»X  ITB024 

XNEG-T  XIB02S 

GO  TO  40  ITB026 

20  NSIGN2*»t  XTB027 

XPOS»X  XT  BO  28 

IPOS»X  IT  E  029 

40  IF  (NTXPB.EQ.2)  GO  TO  SO  ITB030 

SO  IF{  (NIT-1)  .Le.O)  GO  TO  60  ITB031 

NSIGN>NSIGN1*NSIGN2  ITH0J2 

IF  (NSIGN.LE.  0)  GO  TO  BO  II8033 

60  NSIGN1«NSIGN2  ITB034 

NIT*  VII*  1  ITB035 

C  INCBEHENT  TO  FIND  SOLUTION  IHTBBVAL  ITB036 

X*X»SIGN«DX  X3B037 

GO  TO  100  ITB038 

C  INTEBPOLATION  FOB  SOLUTION  ITB039 

80  NTXPB»2  XTB040 

NXI-NIT*1  IT  3  04  1 

XSAVE»X  ITB042 

BATIO-  (XPOS-XNEG)  /  (XPOS-XNEG)  ITB0U3 

X»XNBC*BATIO*(IGIVEN-XNEG)  ITB044 

IF  (ABS  (X— XSAVa)  -  BBBOBX)  90,90,100  IXH04S 

90  NTT?*>3  ITB046 

100  BBIUBN  XT B 04  7 

C.ND  ITS  04  8 


n  n  n  n  n  n  n 


UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  ITMS1 


XBS002 

xasoo3 

subbooixbb  nasi  (asi,pa3so,pa3soi,BcouBTJasixo,asiHi)  iasoou 

las  005 

ruis  sobbooiibz  sbabches  fob  ah  asi  a£mis  iub  exit  pxabb  zasoo6 

PB2SS0BB  COBDXTXOB  IB  COBJUMCTXOB  XXXH  OCV.  BsUs.aE*aB.  XBS0G7 

IBS  008 


BBAi.  asi.asiLO.asiaz 

x as 009 

IFtBCOUBT.GT.  11  GO  XO  20 

XBS010 

asi*.  500 

ISS011 

as  ixo-o.  o 

xasou 

as  lax* i.oo 

ISS013 

BC0UBT*2 

xaso  14 

BBXOBB 

xasois 

XBCBBASXBG  BS 1  OBCB EASES  PB3S0  SO: 

xaso  is 

20  XF  (PB330.GT.  PBSSOX)  GO  to  40 

xasoi7 

asiax»asi 

xasois 

GO  TO  60 

XBS019 

40  asixo*asi 

xas  02  o 

60  asi*  (asiux*asixo)/2. 0 

I8S021 

NCOOBT*BCO(1BT  +  1 

XBS022 

BEtOBB 

xas 023 

aBD 

18  3024 

nnn 
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UNSTEADY  CONSTANT  AREA  EJECTOR  PROGRAM  (UCAE) 
SUBROUTINE  SUMUP 


suaOOx 

SUBU03 

3UBB0UTIBB  SUaDPtfP»ESS.BCYC.PSOPOI,CDBtE)  SUBJ04 

C  TILLS  SUBBOUTIBE  Cii.CDi.iZS5  THE  B*iB  tiLUES  OP  SSI.  PSlPI,  B Si,  SUBOOS 

c  paoso,  psopo,  ino  aa3.  it  i  Dili  piss  is  10  be  seed  xi  suould  suaooo 

C  COBTilB  iH  IBTEGEB  BUSSES  OP  PBSSSDBS  CYCLES.  SUB007 

c  suaooa 

coaaoB/BiocK2/BSi.aa3, psipi, paoso. bsp.psjso  suaoo9 

coaaoa/BLOcsJ/sasi.spaoso,saa3,spsopo,spsiPi,SHSP  soaoio 

BESS  BS1.BB3  SUfldli 

LOGIC iX  COBfS  SUB012 

IP (IPBSSS. GT. 1 )  GO  10  10  SUa013 

sas i3  o. o  suaoiu 

5Paoso»o.o  suaois 

saaj-o.o  iuaoib 

spsopo-o.  o  suaoi7 

spsipi»o.o  suaoia 

sksp«o.o  suaoi9 

10  IP  (BCYC.EQ.IPSSSS.  iSD. .SOT .CDBVE)  GO  TO  20  SUB020 

sasi«asUsasi  suboei 

saa3» a  b3»s  bb3  subo22 

S?SOPQ»PSGPOI*SPSOPO  SUB  0  23 

SPS1P1  =  PS1P1*SPS1P  1  SUB024 

SPBOSO-PBOSO+SPBOSO  5UB02S 

SBSP»BSP*SBSP  SUB  026 

IF  (IP BESS.  BE.  BC  PC)  BSXOBB  5US027 

20  TQT=PL01T (BCXC)  5QB028 

sasi*sasi/Toi  sua029 

saa3=saa3/Toi  suao3o 

SPSOPO-SPSOPO/TOT  SUS031 

SPS1P  1*SPS1 PI/TOT  5U8032 

SPaOSO«SPBOSO/TOI  sua 033 

SBSP=*S BSP/TOT  SUB034 

MBITS  {6. 30)  saSI.SSBS.SPSGPO.SPSIPl.SPBOSO.SBSP  SUB035 

30  P0BB1I  (///,10X."TUE  BSiB  ViLUSS  OVEB  i  PBSSSDBS  CYCLE  ifiE:",//,  SUB036 

& 10E,»aS1»».P10.4,10i,«aa3-“, P9.4,11I,MPS0P0-«,r7.4.//, 10X,  SUB037 

i“PSlP1»«,P8-«.  10X,"PaOSO»«,P8.4.  10X ,  "BSP»",  P8.  4 1  S08038 

BSTQBB  S0B039 

SBD  SUB040 


m 


A.  2  SAMPLE  PROGRAM  INPUT  AND  OUTPUT 


-HE  INPUT  IS: 

JIOCAE  HP1  =  2.0,AP1N  3=.3  3  3,  H9AI=3.  3, PH3SOI=2. 34* 

THE  OUTPUT  IS: 

INPOT  DATA  FOB  CASE  NO.  1 

3  S*  1.4000  GP*  1.4000  HNSP*  1.0000 

■"S9P0*  1.0000  AP1H3*  .3330  HP1=  2.0000 

Pfl3SOI=  2.3400 

pOH  TBIS  CASE,  BCRIT*  .8014  PCHIT»  .1951 


P1R  THE 

A  HP* 

SINUSOIDAL 

100 

PRESSURE 

PAYG* 

RATIO 
.  1000 

CURVE: 

NCYC*  10 

TIHE  SLOGS/CYCLE 

SLOG 

TIBS  SLOG  DATA 
BEGINE  PSOPO 

WSP 

ns  i 

PS1P1 

HS3 

PHOSO 

PS3SO 

1 

SB 

.  1059 

.  2*65 

.44  87 

.7215 

.54  35 

2.9384 

2.4036 

2 

HI 

.1095 

.2577 

.4551 

.7434 

.5446 

2.8627 

2.  3398 

3 

HR 

.  1095 

.2577 

.4551 

.7434 

.5446 

2. 8627 

2.3398 

4 

SB 

.  1059 

.2*65 

.44  87 

.7215 

.5435 

2.9364 

2.4036 

5 

SB 

.  1000 

.2205 

.4186 

.6935 

.5362 

3.0744 

2.5280 

6 

SB 

.  0941 

.  1946 

.  3867 

.6643 

.5287 

3.2275 

2.66  79 

7 

SB 

.  0905 

.  1786 

.  3659 

.6454 

.5240 

3.3328 

2.7640 

8 

SR 

.0905 

.1786 

.3659 

.6454 

.5240 

3.3328 

2.7640 

9 

SR 

.0941 

.1946 

.  3867 

.6643 

.5287 

3.2275 

2.6679 

10 

SB 

.  1000 

.  2205 

.4186 

.6935 

.5362 

3.0744 

2.5280 

THt  NEAR  VALOF.S  OV 

EB  A  PRESSORE  CYCLE  ABE: 

NS1* 

.4150 

BH3 

* 

5354 

PSOPO 

«  .1000 

PS  IP  1* 

.6936 

PHOSO*  3 

.0872 

WSP* 

.2196 

THE  TOTAL  NONE EF  0*  TINE  SLOGS  BEAD  NAS  10 
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APPENDIX  B 

The  initial  phases  of  this  research  effort  were  concerned  with  the 
development  of  a  large-scale,  cold-flow,  air  fluidic  oscillator  to  pro¬ 
vide  a  periodic  ejector  driver  flow  and  an  investigation  of  the  potential 
use  of  fluidic  devices  to  enhance  mixing  between  multiple  supersonic 
streams . 

B.l  WATER  TABLE  EXPERIMENTS 

The  preliminary  stages  of  this  investigation  were  based  on  the 
hydraulic  analogy  and  were  conducted  entirely  on  the  water  table.  With 
this  approach,  simple  and  inexpensive  experimental  models  could  be  rapidly 
fabricated  which  would  provide  a  qualitative  insight  into  the  factors 
influencing  the  flow. 

The  use  of  the  hydraulic  analogy  is  generally  accepted  for  modeling 
one-dimensional  or  plane  two-dimensional,  compressible  flow;  however,  the 
analogy  for  unsteady,  compressible  flow  modeling  is  not  as  well  understood 
or  accepted.  Loh  [40]  derives  the  constitutive  equations  governing  one¬ 
dimensional  channel  flow  of  a  liquid  and  makes  a  comparison  with  the 
corresponding  equations  governing  unsteady,  compressible,  isentropic  flow. 
He  concluded  that  qualitatively  the  analogy  exists  although  the  relation¬ 
ships  between  the  two  flows  are  considerably  more  complex. 

All  preliminary  water  table  experimentation  was  done  with  styrofoam 
models.  These  models  were  cut  with  an  electrically  heated  nichrome  wire 
from  four  inch  thick  styrofoam  blocks.  To  improve  motion  pictures  of  the 
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flow,  a  black  polymer  film  was  bonded  to  the  styrofoam  blocks  prior  to 
cutting  so  that  the  upper  surface  of  all  models  would  be  highly  visible 
against  the  white  bedplate  of  the  water  table.  The  styrofoam  models  were 
mounted  on  a  wooden  bedplate  with  specially  fabricated  spiked  clamps. 

The  bedplate  was  clamped  in  position  on  the  water  table  surface;  gates 
and  ports  on  the  bedplate  controlled  the  flow  of  water  to  the  oscillator 
models. 

Three  basic  fluidic  oscillator  geometries  were  investigated  and  are 
shown  in  Figs.  B.l-1,  B.l-2,  and  B.l-3.  The  principles  of  operation  of 
the  single-jet  oscillator  are  discussed  fully  in  Section  3.1;  in  a  similar 
fashion,  the  balance-ended  oscillator  exploits  the  bistability  of  a  jet 
impinging  on  a  knife-edge.  The  middle  jet  of  the  three-jet  interleaved 
oscillator  is  steady  and  is  supplied  through  a  port  on  the  surface  of  the 
water  table  bedplate.  These  oscillator  configurations  were  selected 
because  of  their  simplicity  and  adaptability  to  the  problems  of  interest. 

A  flip-flop  jet  oscillator  similar  to  Viets1  configuration  [30,31]  was  not 
considered  because  of  the  potential  for  high  flow  losses  and  its  relative 
complexity. 

The  experimental  procedure  during  the  hydraulic  analogy  experiments 
was  to  gradually  increase  the  height  of  the  water  in  the  large  plenum 
upstream  of  the  oscillator  model  being  investigated  until  a  desired  level 
was  reached  and  then  the  observations  were  made.  The  individual  styrofoam 
pieces  of  the  model  could  be  repositioned  easily  so  that  the  effects  of 
the  relative  positions  of  the  pieces  could  be  investigated  quickly. 
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The  objectives  of  these  preliminary  hydraulic  analogy  experiments 

were: 

1.  To  investigate  the  operational  characteri sties  of  these 
oscillator  configurations; 

2.  To  determine  the  factors  influencing  their  behavior; 

3.  To  determine  if  the  pulsations  of  the  three-jet  oscillator 
enhanced  mixing  between  interleaved  supersonic  streams; 
and 

4.  To  examine  a  supersonic-subsonic  ejector  configuration  on 
the  water  table  which  is  driven  by  the  single-jet 
oscillator  in  order  to  gain  insight  into  the  flow 
phenomena  involved. 

The  single- jet  oscillator  configuration  shown  in  Fig.  B.l-1  performed 
very  poorly.  The  motion  of  the  jet  was  extremely  irregular,  moving  only 
slightly  left  or  right  of  the  knife-edge  and  producing  nearly  negligible 
pressure  pulses.  The  centerbody  forming  the  feedback  loop  was  recontoured 
as  shown  in  Fig.  B.l-4.  The  tendency  of  a  jet  to  attach  to  a  nearby  wall, 
the  Coanda  effect,  tends  to  pull  the  inlet  jet  away  from  the  knife-edge 
so  that  it  is  directed  into  the  feedback  loop.  The  narrower  feedback  loop 
channel  and  the  now  vigorous  action  of  the  jet  resulted  in  a  strong  pres¬ 
sure  pulse  moving  through  the  feedback  loop.  The  proximity  of  the  feed¬ 
back  loop  centerbody  to  the  inlet  jet  also  causes  the  inlet  jet  to  flip 
well  across  the  knife-edge  into  the  resonance  chamber  when  the  feedback 
loop  pulse  reaches  the  inlet.  This  is  because  the  surge  of  flow  in  the 
feedback  loop  can  exit  the  loop  only  by  deflecting  the  inlet  jet  toward 
the  resonance  chamber. 
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Figure  B.l-5  is  a  photograph  of  a  plexiglass  model  with  the  feedback 
loop  centerbody  moved  away  from  the  inlet  jet.  Figure  B.l-6  is  the  same 
model  with  the  centerbody  moved  slightly  towards  the  jet.  The  jet  action 
in  the  latter  figure  is  much  more  pronounced  and  stable. 

The  balance-ended  oscillator  exhibited  stability  problems.  If  the 
stagnation  height  of  the  water  in  the  plenum  was  gradually  and  carefully 
increased  during  start-up,  the  oscillator  operated  at  about  .44  Hz.  A 
rapid  start-up  of  the  oscillator  or  a  disturbance  in  the  oscillator  would 
result  in  the  oscillator  operating  at  approximately  twice  the  previous 
frequency,  about  .90  Hz,  and  at  half  of  the  previous  amplitude.  Apparently 
the  lack  of  surfaces  near  the  jet  leave  it  free  to  oscillate  in  higher 
frequency  modes.  If  the  resonance  chambers  are  capable  of  sustaining 
these  higher  frequency  modes,  the  oscillator  will  be  somewhat  unpredictable 
and  unstable  operation  will  be  observed.  Control  vanes  were  positioned  in 
the  vicinity  of  the  jet  inlet  as  shown  in  Fig.  B.l-7.  The  tendency  for 
the  jet  to  attach  to  these  surfaces  caused  the  jet  to  flip  further  to  the 
left  and  right.  This  action  prevented  the  jet  from  oscillating  in  higher 
frequency  modes  and  increased  the  amplitude  of  the  pulses. 

Figure  B.l-8  is  the  flow  field  produced  by  the  three-jet  oscillator. 
The  large  undulations  in  the  steady  jet  at  the  center  were  absent  when  all 
three  jets  were  steady. 

The  single-jet  oscillator  was  also  used  to  drive  a  supersonic-subsonic 


configuration  during  the  hydraulic  analogy  experiments.  Not  much  time  was 
spent  on  this  model  due  to  leakage  and  secondary  flow  control  problems. 
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The  effects  of  the  pulsed  driver  on  the  shock  wave  structure  within  the 
mixing  duct  and  the  rapid  mixing  of  primary  and  secondary  streams  after 
the  disintegration  of  the  supersonic  core,  see  Section  3.3.2,  were  observed. 

A  Super-8  color  film  detailing  the  above  experiments  is  available 
upon  request  from  Professor  A.  L.  Addy. 

The  final  series  of  hydraulic  analogy  experiments  were  performed 
with  a  plexiglass  model  of  the  single-jet  oscillator;  this  is  shown  in 
Fig.  B . 1 -9 . f  The  plexiglass  model  has  greater  dimensional  stability  and 
control  than  the  styrofoam  model.  The  experiments  with  the  styrofoam 
model  gave  only  a  crude  picture  of  what  factors  affected  oscillator  per¬ 
formance  and  how.  Features  incorporated  into  the  plexiglass  model  included 
(1)  adjustable  inlet  and  outlet  jet  widths,  (2)  adjustable  inlet  to  knife- 
edge  separation  distance,  and  (3)  a  movable  feedback  loop  centerbody. 
Leakage,  which  had  been  a  problem  with  the  styrofoam  model,  was  prevented 
by  using  o-rings  and  a  vacuum  hold-down  system. 

Frequency,  waveform,  and  amplitude  characteristics  of  the  single- jet 
oscillator  model  were  accurately  determined  by  videotaping,  through  the 
transparent  glass  and  plexiglass  walls  of  the  water  table  and  resonance 
chamber  of  the  oscillator  model,  the  passage  of  waves  in  the  resonance 
chamber.  A  ruled  scale  and  timer  were  located  within  the  camera's  view. 

Stop  action,  single-frame  advance,  and  slow  motion  features  of  the  video 
system  aided  accurate  waveform  and  amplitude  determination. 

fThis  series  of  experiments  and  the  design  of  the  plexiglass  oscillator 
model  were  performed  by  T.  R.  Stover  [42],  former  graduate  student. 
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A  series  of  experiments  were  conducted  with  this  oscillator  model 
following  a  classical  experiment  sequence.  The  variables  of  interest  were 
the  ratio  of  the  outlet  jet  to  inlet  jet  width,  the  ratio  of  the  inlet  to 
knife-edge  separation  distance  to  the  inlet  jet  width,  and  the  ratio  of 
the  inlet  jet  width  to  the  oscillator  resonance  chamber  length.  These 
ratios  will  be  referred  to  as  R} ,  R2 ,  and  R3 ,  respectively. 

Oscillator  frequencies  were  found  to  be  12%  less  than  what  surface 
wave  theory  suggests.  This  is  certainly  due  in  part  to  the  time  required 
for  the  jet  to  flip  across  the  knife-edge  and  then  back  again.  Also,  the 
ratio  of  the  wave  amplitude  to  mean  water  height  was  not  negligible. 

Amplitudes  were  found  to  rapidly  increase  with  the  ratio  of  knife- 
edge  separation  distance  to  inlet  jet  width,  R2 ,  to  some  maximum  value 
and  then  to  gradually  decrease  as  this  ratio  was  further  increased.  The 
value  of  R2  at  which  the  maximum  amplitudes  occurred  was  insensitive  to  Rt 
but  increased  slightly  as  R3  was  decreased. 

The  ratio  of  the  mean  height  of  the  water  in  the  resonance  chamber  to 
the  plenum  chamber  stagnation  height  was  linearly  proportional  to  Rt .  Wave 
amplitudes  were  found  to  increase  with  Rj  up  to  Rj  =  1.0.  Lesser  ampli¬ 
tudes  were  always  observed  when  Rx  was  further  increased  to  Rj  =  1.3. 
Generally  amplitudes  were  found  to  increase  with  R3  too.  These  facts  sug¬ 
gest  that  the  maximum  amplitudes  are  obtained  at  the  maximum  mass  flowrate; 
i.e.,  when  the  inlet  jet  has  the  most  energy  as  it  flip-flops.  Choking  at 
the  inlet  appears  to  occur  if  Rt  >  1.0. 


The  waveform  was  not  significantly  altered  by  changes  of  any  of  these 
geometric  ratios.  The  waveform  rise  times  decreased  slightly  as  the  mass 


flowrate  was  increased.  A  higher  frequency  beat  was  always  noticed  at  the 
peak  of  each  pulse.  This  was  due  to  transverse  oscillations  of  the  jet 
while  it  resided  in  the  resonance  chamber.  This  is  believed  to  be  the 
cause  of  the  high  frequency  beats  seen  in  the  air  flow  oscillator,  see 
Fig.  3.3.2-29. 

The  results  of  these  experiments  were  used  to  design  the  air  flow 
oscillator.  A  value  of  Rx  was  selected  that  was  sufficiently  small  to 
ensure,  beyond  any  doubt,  that  the  C-D  primary  flow  nozzle  at  the  ejector 
inlet  was  always  choked  but  Rj  was  still  sufficiently  close  to  a  value  of 
1.0  such  that  the  maximum  amplitude  pressure  pulses  were  approached; 

Rt  =  2/3  was  used.  The  oscillator  design  incorporated  a  movable  knife- 
edge  so  that  amplitudes  could  be  adjusted.  It  was  assumed  that  the  mass 
flowrate  through  the  oscillator  would  be  large  enough  so  that  the  optimum 
value  of  R2  would  not  be  too  large.  A  moderate  value  of  R2  was  used  to 
determine  the  maximum  knife-edge  separation  given  the  inlet  jet  width. 
Unless  this  value  of  R2  was  far  too  small,  the  amplitude  of  the  pulses 
obtained  at  the  maximum  knife-edge  separation  would  be  near  the  maximum 
obtainable  for  the  previously  selected  value  of  Rt  . 

Details  concerning  the  design  and  performance  of  the  large-scale  air 
flow  oscillator  can  be  found  in  Sections  3.1  and  3.3  of  this  report. 


B.2  AIR  FLOW  FLUIDIC  OSCILLATOR  EXPERIMENTS 

A  series  of  cold  air  flow  experiments  were  run  with  the  fluidic 
oscillator  prior  to  all  ejector  experiments.  Frequencies  from  75  Hz  to 
250  Hz  were  examined;  the  highest  frequency  was  limited  by  the  minimum 
possible  lengths  of  the  feedback  loop  and  resonance  chamber. 
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The  operating  frequencies  of  the  oscillator  were  about  25%  less  than 
the  estimates  made  assuming  pulses  moved  at  the  estimated  sonic  speed. 

By  changing  the  resonance  chamber  length  while  keeping  the  feedback  loop 
length  constant  and  measuring  the  resulting  frequencies,  it  was  possible 
to  estimate  pulse  propagation  speeds.  Propagation  speeds  were  in  excess 
of  the  sonic  velocity  as  expected.  Thus,  it  appears  the  time  it  takes  the 
inlet  jet  to  flip-flop  is  the  cause  for  the  low  operating  frequency  limits. 

The  ratio  of  the  resonance  chamber  to  feedback  loop  lengths  were 
varied  over  a  wide  range  with  no  apparent  effects  other  than  frequency 
effects. 

The  results  of  these  experiments  were  used  to  determine  which  operat¬ 
ing  frequencies  and  oscillator  configurations  would  be  used  in  the  ejector 
experiments.  The  knife-edge  setting  was  left  at  the  maximum  R2  or 
separation  distance.  Pressure  pulses  were  considered  to  be  sufficiently 
large  and  extending  the  knife  could  have  caused  vibration  or  leakage 
probl ems. 
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Figure  B.l-4  The  improved  single-jet  oscillator  configuration 
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ure  8.1-5  Single-jet  oscillator  with  centerbody  moved  away 
from  inlet  jet 
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Figure  B.l-6  Single-jet  oscillator  with  centerbody  in  position 
for  jet  attachment 


Figure  B.l-8  Three-jet  interleaved  oscillator  flow  field 
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